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SECTION 1 


INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

■k 

This report is the Third Interim Technical Report covering three 
studies performed by The Charles Stark Draper Laboratory, Inc. , for 
the NASA Goddard Space Plight Center during the period from August 1, 
1974 to October 31, 1975. These studies were as follows: 

(1) Modified Landmark Utilization Study - An investigation into 
the use of range and range-rate measurements combined with 
known landmark sightings using earth sensor imagery to esti- 
mate spacecra-^t attitude, orbital e’jhemeris, and gyro bias 

f 11 ** 

drift. The original Landmark Utilization Study had 

investigated the combined use of star and known landmark 
sightings . 

( 2 ) Development of Star and Landmark Measurement Equations - 

A detailed analysis of star and landmark measurement errors 
and sensitivities. 

(3) Horizon Utilization Study - An investigation into the use 
of horizon measurements for estimation of spacecraft atti- 
tude and gyro bias drift. 

Most of these studies represent a natural follow-on to previous 
studies performed by CSDL for GSFC pertaining to the optimal treatment 
of attitude and orbital ephemeris information for spacecraft that gen- 
erate high-resolution imagery of the Earth. (A brief description of 
the previous studies is given in the next subsection) • One of the long- 
range goals of this overall effort has been to determine how spacecraft 
attitude and orbit data can best be used to improve the mapping accur- 
acy of a multi-spectral scanner. 

* A fourth task involving the actual use of array searchlights with 
Landsats 1 and 2 is still underway and a complete report on this 
effort will be made in the Pinal Report. 

** Superscripts refer to numbered references at the end of this report. 
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1.2 PREVIOUS STUDIES 


Since July 1971, the Draper Laboratory has performed various 
studies for GSFC pertaining to the determination of attitude and orbital 
ephemeris for an advanced Earth Observation Satellite (EOS) . These 
are briefly described chronologically in the following paragraphs. 

Initial SIMS Trade Study . This study provided GSFC 

with technical data thar could be used for selection of an optimum 
Stellar-Inertial Measurement System (SIMS) for EOS. From a large 
number of initial SIMS candidates the following three candidates were 
selected for detaixed evaluation; 

SIMS-A Strapdown gyros and 

strapdown star mapper 

SIMS-B Strapdown gyros and 

gimbaled star tracker 

SIMS-C 3-Axis gimbaled gyro 

platform and strapdown 
star mapper 


Derived from 
Honeywell SPARS 


Derived from TRW 
PPCS/PADS 


Subsystems 
defined by CSDL 


A significant accomplishment of this study was the generation of 
statistical data showing how well the spacecraft attitude and gyro bias 
drift could be estimated for each candidate. Use was made of an opti- 
mal smoother to obtain 'least squares' estimates of attitude and gyro 
bias drift after processing star measurements over several orbits. 

( 5 ) 

Detailed SIMS-A Study . This study was a more detailed and 
accurate Monte Carlo simulation of SIMS A. This confirmed the statis- 
tical estimates of accuracy obtained in the previous study. Performance 
data v/as also generated showing the effects of certain bias errors that 
had not been modeled before in earlier covariance studies. These in- 
cluded errors such as gyro scale factor error and gyro input-axis mis- 
alignment. 
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Attitude Determination with Landmarks . This was an investigation 

determine how well spacecraft attitude and gyro bias drift could 
be estimated using known landmarks in the sensor imagery. The pro- 
mising results indicated a possible backup method of attitude deter- 
mination in the event of star sensor failure. 

Mapping Error Sensitivity to Errors in Spacecraft Attitude and 
Orbital Ephemeris . The main objective of this study was to deter- 
mine the effects of errors in spacecraft attitude and in orbital ephe- 
meris on the ability to determine the locations of unknown landmarks 
relative to known landmarks. Emphasis was placed on landmark location 
accuracy within the continental USA. One significant result was the 
discovery that errors in spacecraft attitude and gyro bias drift have 
a reduced effect on mapping errors due to the presence of certain 
negative correlations. 

Attitude and Ephemeris Determination with Landmarks and Stars . 

This study was an investigation into the use of star sightings 
together with known landmarks to estimate spacecraft attitude, gyro 
bias drift and orbital ephemeris. A 12-state covariance analysis was 
employed for the simulations. One important finding v/as that known 
landmarks could be used together with star sightings to satisfactorily 
estimate orbital ephemeris. 

Artificial Landmark Implementation . This study was an investi- 
gation into the feasibility and practicality of establishing a national 
system of artificial landmarks suitable for automated recognition. 

I . 3 SUMMARY 

Section 2 reports the results of the modified landmark utilization 
studies. A covariance analysis is used to compute satellite state 
uncertainties with emphasis on systems involving xandmark sightings, 
and tracking facilities. Equations governing the coordinate trans- 
formations, the stats dynamics, and the estimation process are given. 
Simulation results are given for parametric variations of measurement 
noises and update frequencies, and for various combinations of measure- 
ment systems. 
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Section 3 derives generalized measurement equations for star and 
landmark sightings. The derived equations are specific to certain types 
of star trackers and landmark sensors, but consider general types of 
noise sources. 

Section 4 presents the results of a stu<ly aimed at the use of 
infrared horizon measurements for satellite attitude and gyro bias 
drift updates. Measurement equations are derived, the state and esti- 
mation equations are given, and simulation results are presented. The 
results consider variations in measurement noise and update frequency 
for both sun-synchronous and earth-synchronous orbits. 
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SECTION 2 


MODIFIED LANDMAEK UTILIZATION STUDY 


2,1 INTRODUCTION 

The purpose of the Modified Landmark Utilization Study is to 
investigate the use of ground or navigation satellite tracking measure- 
ments of range and range-rate together with the use of known landmarks 
observed in the imagery of a multi-spectral scanner (MSS) , to estimate 
spacecraft attitude, orbital ephemeris and gyro bias drift. This sys- 
tem configuration was proposed as an alternative to the one studied 
in the original Landmark Utilization Study The original study con- 
sidered the use of star measurements together with known landmark 
observations. The present study is thus the original study "modified" 
to incorporate range and range-rate measurements. 

This study considers radar measurements of the range and range- 
rate between an Earth Observation Satellite (EOS) and either a ground 
tracking station or a navigation satellite of the TDRSS (Tracking and 
Data Relay Satellite System) type. For the ground tracking configu- 
ration, the use of two unified S-Band System radars located at Merritt 
Island, Florida and at Goldstone, California, are found to be suffici- 
ent to give satisfactory tracking of EOS over the continental USA. For 
the navigation satellite tracking configuration, two geosynchronous 
TDRSS satellites located over the equator at 41 and 171 degrees west 
longitude are employed. 

The Earth Observation Satellite is assumed to be in a circular 
sun-synchronous orbit with an inclination of 99 degrees and an altitude 
of 1000 km (540 nmi) . This study also includes the effect of gravitat- 
ional harmonic uncertainties. Figure 2-1 illustrates a typical obser- 
vational pass over the continental USA. The ground tracking stations 
at Merritt Island and Goldstone are indicated. Shown on board the 
satellite is a multispectral scanner vjhose beam sweeps back and forth 

across the ground track to generate a swath of imagery 90 nmi (167 km) 

* 

wide. 


*Although the svath width for Landsat and the planned EOS is 100 nmi 
(185 km) , the present study is restricted to using landmarks over 
a smaller swath width since this was the desire of GSFC . 
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Under nominal attitude conditions (that is, with the X^-axis 
along the local horizontal, and zero roll and yaw angles) , the space- 
craft body axes (Xg,Yg,2g) are oriented as shown in Figure 2-1, with 
Xg in the direction of spacecraft motion and pointing to nadir. It 
is assumed that the spacecraft has three orthogonal body-mounted gyros 
aligned with its principal body axes. 

In this study optimal filtering techniques are used to process 
the landmark and the range and range-rate measurements. Since the 
study is primarily concerned with obtaining a statistical measure of 
performance in estimating spacecraft attitude, orbital ephemeris, and 
gyro bias drift, a covariance type of analysis is used. Here, the 
covariance matrix of the uncertainties in the estimates of the state 
parameters is propagated and updated using Kalman filter equations. 

This approach yields a statistical indication of achievable performance 
for a given system error and measurement noise model. 

It should be noted that there are drawbacks to limiting an error 
analysis to covariance matrix processing alone. One drawback, of 
course, lies in the simplifying assumptions relative to modeling of 
spacecraft orbit and attitude history. Another is found in the incom- 
plete modeling of system biases and random processes. To minimize the 
effects of modeling inadequacies, all important error sources should 
be modeled. If the statistics for certain states are truly Gaussian, 
the Kalman filter in linear form can be shown to be optimal. However, 
if the statistics of some states are poorly known or incorporate non- 
zero biases, the effort to incorporate poorly modeled states may result 
in unacceptable performance. 

Consider States 


One example of a state modeling problem in the present study is 
that associated with the TDRSS satellite ephemeris uncertainties. Dur- 
ing the relatively short time of data processing in the present study 
(typically, ?.5 2 -. 0 'irs), it would seem reasonable to expect very little 
change in the TDRSS satellite ephemeris uncertainties. Hence, 
measurements by TDRSS can be considered 
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highly time-correlated or as ha\^ing bias components. Rather than 
model these uncertainties as white noises, or estimable biases , another 
approach is to model them as "consider" states. This approach enables 
the Kaiman filter to recognize the presence of the "consider" state 
variable errors which are unobservable while minimizing the number of 
estimator or covariance states. The consider state equations are given 
in Section 2.7. Results are given in Section 2.8.4 comparing esti- 
mation uncertainties with TDRSS ephemeris uncertainties modeled as 
white noises and as consider states. 

Gravitational Harmonic Uncertainties 

The effects of gravitational harmonic uncertainties are included 
in all covariance runs made for this study. Also taken into account 
is a 0.2 PPM uncertainty in the main gravitational constant, u. The 
term-by-term difference between the two most recent Smithsonian gravi- 
tational models, SEIII (1973) and SEII (1969), as represented in the 
C and S expansion coefficients (through n,m =10, 10 ) is initially 
defined to be equal to the one sigma uncertainty in knowledge of these 
coefficients. These differences are subsequently scaled down to match 
propagation data obtained from GSFC for 5 orbits. The scaling effect- 
ively defines the one sigma uncertainties as one-third of the difference 
between SEII and SEIII coefficients. 

To model these uncertainties directly using Monte Carlo techni- 
ques would be excessivly time-consuming. As an approximace and fea- 
sible alternative, ephemeris uncertainty covariance matrices, repre- 
senting the portion of the orbit that is inaccessible to radar track- 
ing, are generated by Monte Carlo runs with randomized gravitational 
constant uncertainty. In the covariance runs made for this study this 
ephemeris uncertainty covariance matrix is incorporated into the simu- 
lations at the start of each pass over the USA, Thus the gravitational 
harmonic uncertainties acting over that portion of the orbit which is 
beyond radar coverage is approximately accounted for. The effect of 
these uncertainties on that portion of the pass over the USA having 
radar coverage is not modeled, since it is relatively negligible. 



2 . 2 SUMMARY 


This report presents the analytical equations and performance 
results of the Modified Landmark Utilization Study. The study investi- 
gates the use of range and range rate measurements by ground or TDRSS 
tracking stations , in combination with landmark and/or star measure- 
ments , to update a state vector which includes EOS spacecraft attitude, 
oribtal ephemeris and gyro bias drift. Optimal Kalman filter techniques 
were used to process the different types of measurements. A covar- 
iance analysis that gives a statistical indication of performance in 
estimating spacecraft attitude and orbital ephemeris was adopted for 
this study. For studies involving TDRSS ephemeris uncertainties, 
"consider" states were added to the covariance data processing to im- 
prove the modeling- The study is restricted to the use of landmarks 
and ground tracking stations within the continental USA and Alaska. 

The effect of gravitational harmonic uncertainties is also incorp- 
orated into the simulations. 

Comparative data is presented showing the performance obtained 
with various combinations of the indicated measurements. Data is also 
presented showing the sensitivity of performance to variations in the 
following parameters: 

(1) Initial uncertainties in spacecraft attitude, orbital 
ephemeris and gyro bias drift. 

(2) Ground tracking station location uncertainty. 

(3) Range and range-rate measurement noise. 

(45 Number of range and range-rate updates per pass. 

(5) Number of landmark updates per pass. 


(6) Landmark position uncertainty. 


The data presented shows that range and range-rate measurements 
alone can be used to satisfactorily estimate orbital ephemeris, but 
that landmark measurements must also be incorporated to satisfactorily 
estimate spacecraft attitude and gyro bias drift. 

2.3 COORDINATE SYSTEMS 

This section defines the reference coordinate systems or frames 
used in the error study and simulations. These primary reference 
frames are the following: 

Basic Inertial (I-frame) 

Orbit-Oriented Inertial (O-fraroe) 

Body-Fixed (B-frame) 

2.3.1 Basic Inertial Coordinate System (I-Frame ) 

The coordinate axes for this system are defined in Figure 2-2. 

The axes and both lie in the equatorial plane with pointing 
tovjards the vernal equinox. Axis 3^ points along the north polar axis 
of the Earth. Star catalogs noimially give the directions of stars in 
this coordinate system. 

2.3.2 Orbit-Oriented Inertial Coordinate System (0-Frame) 

This system of axes is also defined in Figure 2-2. The coordinate 
system is oriented relative to the basic inertial coordinate system 
through the angles n and i. The first angle is the right ascension of 
the orbit ascending node, and the angle i is the orbit inclination. 

This orbital plane does precess slowly about the earth's rotational 
pole due to oblateness of the earth. However, the orbit-oriented co- 
ordinate system is defined herein to be an inertial frame since, in our 
simulations, orbit plane rotation due to precession is ignored as it 
would add unnecessary complexity to the simulations. Since real rota- 
tion of the EOS orbit plane is a small fraction of a degree over the 
course of a typical simulation, the distribution of available stars is 
not affected by such precession. The transformation matrix from 
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(Vernal 

Equinoj^ 






basic inertial to orbit-oriented inertial coordinates is given by: 


— 


— 


— 



1 

0 

0 


cfJ 

sfi 

0 

0 

ci 

si 


-sfi 

cfi 

0 

0 

-si 

ci 


0 

0 

1 


( 2 - 1 ) 


where c denotes cosine , and s denotes sine. Thus a star vector s^ in 
the orbit-oriented frame can be computed, given the star vector in 
basic inertial coordinates. 




= T ^ 
01 -I 


(2-2) 


2.3.3 Spacecraft Body-Fixed Coordinate System (B-Frame) 

The axes of this system are such that X-, Y„, and Z are respect- 

O D O 

ively the roj.1, pitch and yaw axes of the spacecraft. The nominal 
orientation of these axes is as follows: 


X - is along the projection of the spacecraft velocity 
B 

vector onto the local horizontal plane 


Yjj “ is normal to the orbital plane 

Z - is along the local nadir 

B 


The orientation of the B-frame with respect to the 0-frame is shown in 
Figure 2-3. The transformation from the 0-frame to the B-frame is 
through the Euler angle sequence of pitch (0) , roll ((f)) , and yaw (i|() 
as shown in Figure 2-3 and expressed by: 


0 


1 



0 


0 


C(f) 0 -s^ 


c0 s8 0 


T 


BO 


0 


-1 


0 -1 

0 0 


0 ctjj sifj 


0 


1 0 


-S0 c0 0 


0 -si|j G(|) 


s t 0 



0 


1 


(2-3) 
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The input axes for the strapdown gyros are ideally colinear 

with the spacecraft body-fixed axes, so that = Y_ and 

g ti g li 



2.4 SYSTEM EQUATIONS 
2.4.1 General Comments 

In the present study the state to be estimated consists of 
twelve parameters which comprise 3 attitude angles , 3 orbital position 
components, 3 orbital velocity components, and the bias drift of each 
of 3 gyros. These parameters can be represented as the components of 
a twelve dimensional state vector x. The dynamic behavior of the sys- 
tem with the above state can be described by the following nonlinear 
differential state equation: 

X = f (x, u, t) (2-4) 

where t is time and u is the noise introduced by gyro random drift, 
gyro quantization, etc. . A priori information about the statistics 
of the noise u is given by the covariance matrix Q where: 

E [ u(t} = Q . 6 (t-T) (2-5) 

At some arbitrary time, t^, a measurement ^ (t ) is taken. This mea- 
surement is related to the state through the measurement equeition; 

^ 2 - 6 ) 

where v (t^) is the measurement noise of a zero mean white gaussian 
process whose covariance matrix R is given by: 

For the special case where Equations 2-4 and 2-6 are linear, the stand- 
ard Kalman Filter would be used to obtain optimal estimates of the 
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state. However, in the present case, these equations are nonlinear 
and must be linearized before using the so-called "Extended Kalman 
Filter". This linearization is generally performed around the esti- 
mated values of the state. In the present covariance analysis, how- 
ever, there is no estimated state and use is made of the known nominal 
state for this purpose. 

For the Extended Kalman Filter, the new state and measurement 
vectors are defined as: 


(Sx 


X - X. 


-nom 


( 2 - 8 ) 


6z - z - z 

— —nom 


The corresponding state and measurement equations are; 



•Sx 

= F 6x 

+ G 

6u 



(2-9) 


•Sz 

= H 6x 

+ A 

6v 



(2-10) 

matrices 

F, G, H, and 

A are : 





F = 

af 


G 

= 

/I ^ 





-nom 



9 ii ' 

1 

1 

—nom 



! 

u = 0 




u = 0 

(2-11) 

H = 

a h 


A 

= 

9 h 




Tx 

—nom 



a“v" 

^nom 




V = 0 




o 

II 

>f 



2.4.2 State Dynamics 

The original nonlinear state of the system x is defined as the 
following twelve dimensional vector: 


2-10 


( 2 - 12 ) 


= (0, <)., By, X, y, z, v^, Vy, v^) 

where 

0,^,4^ - Euler angles defining the body attitude with 
respect to the orbital frame. 

B - Bias drift of roll (X), pitch (Y) , and yaw (Z) gyros 

X y z 

x,y,z - Components of the spacecraft position vector in the 
Basic Inertial frame 


V , V , V 

x' y' z 


- Components of the spacecraft velocity vector 
in the Basic Inertial frame 


The eqi-ations that relate the Euler rates (0f4'/4*)ho the spacecraft 
body rates ((ij^,ojy,u^) , measured by the X, Y, and Z gyros, are given byi 


r . 1 


— 


— 


— — 

e 


SlJj/ctf) 

-c4j/c<J) 

0 


0) 

X 



cij/ 

S(|) 

0 


w 

y 

‘I. 


S(j)Slp/C(i) 

-S(|)c4i/c(j) 

-1 




(2-13) 


or in matrix notation: 


i = G (i) m (2-14) 

where ^ 

T 

(iJ = [ W , 111 , 01 ] 

— x' y' z-* 

G (^) = the matrix in Equation 2-13. 


The relationship between the gyro measured body rates (w) and the 
true body rates (w^. ) sensed by the gyros is: 


w = iii.. t D 


(2-15) 


m 

where 

— t 


( 0 , 0 ) 


(ii)^ = orbital rate) (2-16) 
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and p represents the gyro drift rate which can be expressed as follows: 


D = B + (2-17) 

where B = gyro bias drift 

n = additive white noise (gyro random drift) . 

The gyro random drift is assumed to be a zero mean white gaussian 
noise whose covariance matrix Q is given by: 

S [n(t) n'^(T)l = 0 • 6{t-r) (2-18) 

Equations 2-14 can therefore be rewritten as: 

k = G(€) (-^ + 0(5) B + G(C)n (2-19) 

* 

Since the components of the B vector are unknown constants, they should 
be included as additional elements to be estimated in the state vector. 
This can be done by augmenting Equation 2-19 and including B in the 
state. Recalling that for constants there are no dynamics: 


= ( 0 , 0 , 0 ) 


( 2 - 20 ) 


Hence, Equation 2-19 can be changed to: 


"* - 



“ 





i 


°3 

G(e) 


r 

+ 

0 (V 

B 









where 0^ is a 3 x 3 null matrix. 

Since the orbital dynamics are uncoupled with the attitude angle 
dynamics, their description is given by a separate set of first-order 
ordinary differential equations: 

^^r the purposes of this study, the gyro bias drift is assumed to be 
unkno^m. However, in normal practice, a portion of the gyro bias drift 
is usually known and is compensated for when integrating the outputs 
of the gyros. 
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where 


v’’ - (v^, Vy, v^) 
T 

r = (x, y, z) 


{ 2 - 22 ) 


)i = gravitational constant 


The combined state of the system for attitude, orbital ephemeris, and 
gyro bias drift can be written as: 


= (r 


B 


V ) 


(2-23) 


where the dynamics of this state are given by Equations 2-21 and 
2-22. However, both equations are nonlinear aid must be linearized 
around the nominal state values. Let a new 12 dimensional state be 
defined as follows: 

m fp rp m m 

Sx = (6? , 6B^, 6r 6v ) (2-24) 


where 


« = 1 - Sr = r - 

B-B,fiv=V“V 
— — —o' — — -o 


(2-25) 


In Equations 2-25 the vectors r and v are the position and velocity 

—"O *"■'0 

for a nominal circular orbit with an orbital rate and 


= (o^t, 0, 0) , 


B 

-O 


= ( 0 , 0 , 0 } 


To linearize Equation 2-21 define: 


F = — 

1 B C 


GU) 03, 
^ “t 




0 

-01 

c 

0 


(2-26) 


= G {£) 


5 = £. 


0 

1 

0 


-1 

0 

0 


0 

0 

-1 


(2-27) 


so that the linearized expression for Equation 2-21 can be written as: 




“ 

— 


— - 


— 



^1 

‘^l 


S| 


^1 

* 






+ 


SB 


0, 

0, 


SB 


0-> 



3 

3 




3 


To linearize Equation 2-22 define: 


(2-28) 


F 


2 


3 -U r 



(r^ - 3x^) 

-3xy 

-3xz 

-3xy 

. 2 - 2 1 
(r -3y ) 

-3yz 

-3x z 

-3yz 

(r^-3z 


where x, y, and z are the components of r^ and r is the magnitude of 
r , The linearized expression for Equation 2-22 can now be written as 


— 


_ 


— 

(Sr 


0. I- 


Sr 


— 

3 3 



6v 


F„ 0- 


Sv 

_ _ 


2 3 




where I^isaSxS identity matrix. 
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Having determined the linearized state equations for both atti- 
tude and orbital ephemeris {Equations 2-28 and 2-30) , these equations 
may be combined to represent the state equation of the full linearized 
state as follows: 

6 k = F 6 x + G 6 ri (2-31) 


where 6 x is given in Equation 2-24 and 


““ ' — 


— — 

Fl Gj^ O 3 0 , 


^1 

O 3 O 3 O 3 O 3 



°3 °3 °3 ^3 

tl 
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°3 

O 3 O 3 Fj O 3 


1 

0 

1 


Since and are time invariant, because of linearization 
around the nominal state instead of the estimated state, a transition 
matrix can be derived for the attitude state Equation 2-28 using 
Laplace transform techniques: 




att 


(k,k-l) 


1 0 

0 

t 

1 

1 

1 — 

0 

< 

0 

0 c 

-s 

I 

0 (l-C)/m^ 

0 s 

c 

t 

I 

I 

(l-C)/m^ 0 

- 03 


I 

I 

I 

"3 


(2-33) 


where k and k -1 denote the times (tj^ and of present and 

previous measurements: is the orbital rate; At^^ = tj^-tj^_^ ;S=sin (X) , 

C = cos (X) ; and X = w At, , 

o K 

The derivation of the transition matrix for the orbital ephemeris 
state equation (Equation 2-30) is more complicated since F„ is time 

, . fg 1 

varying. However, an approximate solution ^ ^ in local vertical co- 
ordinates is : 
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(2-C) 

s 

0 

S/a>o 

2(l-C)Ao^ 

0 


(2S-3A) 

(2C-1) 

0 

2{C-1)A.^ 

(4S-3X)/ti>^ 

0 


0 

0 

C 

0 

0 

sAoo 


w (3X-S) 
o 

iUq(1"C) 

0 

(2-C) 

(3X-2S) 

0 

(2-34) 

a)^(C-l} 

-w S 
o 

0 

-S 

(2C-1) 

0 


0 

0 

— m S 
o 

0 

0 

c 



The transition matrix in Equation 2-34 can be expressed in basic in- 
ertial coordinates as follows: 




(2-34a) 


where: 


ILV 


'ILV 


T = T T 
ILV 01 OLV 


■OLV 


C -S 0 

SCO 
0 0 1 


The transition matrix for the combined state for perturbations on atti- 
tude, gyro bias drift, and orbital ephemeris is therefore: 


4i(k, k-1) = 


■^att 




{2-35) 


where Og is the 6x6 null matrix. 
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2.5 MEASUREMENT EQUATIONS 


For this study the measurement data is primarily collected from 
two instruments, the multi-spectral scanner and the tracking radar. 

The first measures the line-of-sight (LOS) to known landmarks while 
the second measures range and range rate from either ground stations 
or navigation satellites. Also included in the study for comparison 

o 

purposes rs measurement data from the star tracker (8 FOV) of the LOS 
to known stars. The measurements are usually assumed to be corrupted 
by white noise. 

2.5.1 Landmark Measurements 

Landmark measurements are obtained by using a multi-spectral 
scanner (MSS) whose beam sv/eeps back and forth across the ground track 
to generate a 90 mile wide swath of ground imagery (see Figure 2-1) . 
Under the nominal attitude conditions considered in this study, the 

O 

maximum angular excursion of the beam from local vertical is 4.8 . The 
direction of the beam in body coordinates is defined by the angles 
and Here is the sweep angle, measured positively about the 

X body axis, between the z body axis and the beam; and is nominally 
zero and is primarily introduced to account for errors in the remaining 
dimension of angular measurement. 

Since the beam is always near the local vertical, the downrange 
and crosstrack errors in landmark position were modeled as equivalent 
angular errors in and 6^^. It should be noted that a full scheme in- 
cluding state estimation would require a more exact error model. The 
above simplified model was considered to be adequate for this covariance 
analysis. No consideration was given to landmark altitude error since 
the system sensitivity to this error is relatively small. 


The unit vector defining the measured LOS to the landmark in 
body coordinates is: 



Th 0 estimated LOS to the landmark in body coordinates is the 
following unit vector: 

= ^BO '^OI 

where r is the estimated spacecraft position vector and 1 is the land- 
mark position vector, both in basic inertial coordinates at the time 
of measurement. 


2-5,2 Range and Range-Rate Measurements 

The range and range-rate measurements represent the range and 
the range rate between EOS and either a ground tracking station (USES 
radar) or a navigational satellite. 

The range from a tracking station to the spacecraft is simply: 

\ ^ " EtI (2-38) 


and the range rate is: 


R_ 


(v - v^) 


(r - 






(2-39) 


where r, v are the spacecraft position and velocity vectors, and 
V.J, are the tracking station position and velocity vectors in basic 
inertial coordinates at the time of measurement. 


For this study the tracking station may be either (1) a ground 
tracking station, or (2) a geosynchronous navigational satellite- For 
either type of tracking station, the station velocity is given by: 





X 


-T 


(2-40) 


where ^ is the inertial earth angular rate of the tracking station 
with respect to the center of the Earth. 
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2.5.3 Star Measurements 


Only one type of star tracker is considered in this study, and 
it is used only to supplement other types of measurements. This is the 

O 

8 FOV star tracker somewhat similar to the CT401 tracker of Ball 
Brothers Research Corporation. This is a strapdown tracker with no 
gimbals. The measured direction of a star passing through the FOV is 
electronically indicated by two angles and 3^,- These are used in 
the following equation to obtain a unit vector defining the measured 
LOS to the star in body coordinates. 



c3^ 

0 

-S8.J, 


1 

0 

0 


0 



% = 

0 

1 

0 


0 




0 

- 

-SUt 


s3t 

0 

c3^ 


0 




-1 


-c3,j,ca,j, 


The estimated LOS to the star in body coordinates is given by 
the following unit vector: 


-B 


= T 


BO 01 




(2-42) 


where s^ is the unit LOS vector to the star in basic inertial coordi 
nates and T„„ is a function of the estimated body attitude angles 6, 
4) and (See Equation 2-3) . A more detailed model of this tracker 
is discussed in Section 3. 


2-5,4 Line arizat i on of Measurement Equations 

The extended Kalman Filter requires that the measurement equations 
be in a linear form, here given as: 


6z = H. 6x + H . 6x + A. (2-43) 

— 1 — Cl — c 1 

where is the sensitivity of the measurement to the state, is 

the sensitivity to the consider states, and is the sensitivity 

of the measurement to the noise. (Consider states are discussed in 
Section 2.6). The subscript i is replaced by M for the MSS landmark 
sightings, by R and R for the range and range rate measurements and by 
T for che star tracker measurements. 
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Landmark Measurement 


Linearization of the estimated landmark measurement equation 

* * 
yields the following; 


'M 




“ 

almg 

fix 

— —0 


t)x 



x=x 
— -o 


_3_ 

3x 


■x=x 
o 


x=x 
— ~o 


(2-44) 


[■fBo] Tqj unit (1-r) + (2-45) 


-e 




0-0 
0 f 0 


) I 

“ I H.. I 

0 ! 


-X , 
t 


H in i » i “ I “ 

-y I 1 0 1 0 I 0 

I ; r t * 


(2-46) 


where the following are the submatrices (or vectors) due to differentia- 
tion with respect to the state elements; 


H. 


-c6 

0 


-5d 

0 


Unit (1-r) 


(2-47) 


0 

-c6 


0 

•s0 


Tqj Unit (1-r) 


(2-48) 


H.r. 


0 

-S0 


0 

C0 


Tqj Unit (1-r) 


(2-49) 


H = 


T T 
BO 2 01 


d^ 


f^x - 

(ly - 

(1, - 


X) 2 -d2 

y) 

3) (Ij^ - X) 


(2-50) 


H = T„- 1 

— y BO 2 01 

d3 


<^x - 
- 

(1„ - 


x) (ly - y) 

y) 2 - d^ 

z) (ly ” y) 


(2-51) 


*Mote that only two rows are given for the matrix since only two 
of the three components of a "unit" LOS vector are needed to complete- 
ly define the line-of-sight to a landmark, star, etc.. 
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1 


(2-52) 


Sz 


T T 

01 

2 d3 


(1^ - x)(l^ - z) 
(ly - y) (1^ - z) 

(1^ - z)2 - 


where 1^^, 1^, and are the components of the landmark position 1 
computed in basic inertial coordinates at the time of measurement: 
TjgQ is a 2 X 3 matrix containing the first two rows of TgQf and 

d=^\l-r\. 


As stated earlier, landmark position errors and MSS measurement 
errors are both modeled as noises in and The sensitivity of the 

landmark measurement to these errors is represented by the following 
matrix : 






M M, 


-ca,,cg,, 
M M 




Meas 




Meas 


(2-53) 


Range Measurement 

Linearization of the range measurement equation yields: 


H 


R 


9x 


<l£ - 


x=x 
o 


x=x 
— . —^1 


= (0, 0, 0, 0, 0,0, 


^^Ry' ^Rz' 


0 , 0 , 0 ) 


(2-54) 


(2-55) 


where the following equations are due to differentiation with respect 
to the state elements: 
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(2-56) 


"Rx 


R 


^^0 " 


Ry 


1 

R 


(ri - rjj) 


Hjjg (r-5 — iTfpj) 


2 T2‘ 


(2-57) 


(2-58) 


Range measurement error is modeled as affected by two types of 
noises^ first by the additive noise on the range measurement itself, 
and second by the noise <^ue to uncertainty in tracking station ephemeris . 


The sensitivity to the first noise is simply unity. 


A 


Rl 


3R. 


9v 


Rl 


= 1 


(2-59) 


The sensitivity to tracking station ephemeris uncertainty modeled 
as white noise is: 


R2 


3R 

^-R2 


9 

9v 


r - r 


£t=^to 


■R2 


(2-60) 


£t“^to 


\2 " %y' ^Rz' 


(2-61) 


where the H’s are given by Equations 2-56 through 2-58. Here the 
change in sign is due to the fact that differentiation is with respect 
to r^ and not to r as in the H equation. 

The corresponding measurement noise variance equation for use in 
the Kalman update equation is ; 


R 


^Rl ^1 *R1 


T 


^R2 ^2 ^R2 


(2-62) 


*Note here that the term "ephemeris" is used for a tracking station 
on the ground or in orbit. 


ft 




J 
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where 


= [Var (RN) 1 
RN 


R_ 


the range measurement noise 
"^IL \ph 


(2-63) 


(2-64) 


and where Rgpj^ is the tracking station ephemeris uncertainty covariance 
matrix in local vertical coordinates and T_t is a transformation matrix 

XJj 

from local vertical to basic inertial coordinates. 

Range Rate Measurement 

When the range rate equation (Equation 2-39) is linearized, the 
result is: 


?R, 


R 


s 

ax 




— 

“i 


- 9 

(V - v^) 

Ifr 

1 

il’ 


l£ - 

Efp f 


(2-65) 


"0 


Hr = 


- ( 0 , 0 , 0 , 0 , 0 , 0 , ( 2 - 66 ) 


where the rollowing equations are obtained by differentiation with 
respect to the state elements: 


^Rx ^2 ^'^0 ^TO ^ ^ 


(2-67) 


^RY “ p2 ^^1 “ ’^Tl “ ^Tl^^ 


( 2 - 68 ) 


" „2 ^^s ^^^2 ^T2 ^“^s ^^2 

^s 


(2-69) 


2-23 


(2-70) 


*RX 


‘^0 “ ^TO> 


‘RY 


R. 


“ ^Tl) 


= 


R. 


'=^2 - >^ T 2 ' 


(2-71) 


(2-72) 


Range rate measurement uncertainty is also modeled as affected 
by two types of noises. First is the noise in the doppler measure- 
ment itself, and second is the noise Vj ^2 due to tracking stacion ephemeris 
uncertainty. 

The sensitivity to the first noise is simply unity. 



3v3 


The sensitivity to tracking station ephemeris uncertainty 


is : 


^2 = 


8R. 


3 \r 
-R2 


3v 


r = r 
-T -To 

V =v 
-T ^0 


•R2 


(v - v^) 






(2-74) 


and 




-T”^T0 




R2 


[H- , H* , H« , H* ■ , H* • , H* • 1 

^ R.x' "Ry' "rz' ”rx' "Ry' 


(2-75) 


where the H’s are given by Equations 2-67 through 2-72, Here the 

change in sign is due to the fact that differentiation is with respect 

to r„ and v and not r and v. 

-T -T - — 
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The corresponding measurement noise variance equation for use in 
the Kalman update equation is: 


” -^Rl ^3 ^R1 ^R2 ^4 *R2 


(2-76) 


where 


= [Var (RRN) ] 


(2-77) 


RRN = the range rate measurement uncertainty 


R* = T-. R . T_^ ■ 
4 IL eph IL 


(2-78) 


where and have been defined previously. 


Star Tracker Measurement 


Linearization of the estimated star tracker measurement equation yields 

A 

the follov.’ing: 


where 


_3 

^21 


BO 


X = Xr 
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1 

1 
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-c9 -s9 0 

0 0 0 


"^01 -I 


IL. 
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0 0 0 
-G0 “S0 0 


"^01 -I 
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0 

-s0 


0 

C0 


1 
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’’^01 -I 
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(2-82) 


The measurement angles of the star tracker are: 


Meas 


+ v 

T rt 

'■True T 


B.. 




+ V, 


Meas 


'True 


*See footnote on page 2-20. 


(2-83) 
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are the random measurement errors which may be ex- 


where v and v„ 


pressed as the components of the random vector 


V 

-ST 


V 




(2-84) 


whose covariance matrix is defined as follows; 

r 2 

0 0 
°^T 

E 6(t - T) 


Tj 


6{t - t) (2-85) 


The sensitivity of the star measurement to the random errors is: 


^-ST 


Ctrr. = a. 


Meas 


sa^ c3^ ca^ 


-ca„ 


= fin 


Meas 
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2.6 KALMAN FILTER EQUATIONS 


The equations for the filterering technique are the extended 
Kalman Bucy equations for discrete measurements. These equations con- 
sist of two parts: 1) propagation between measurements, and 2) updating 
whenever a measurement is taken. Since the linearized state equations 
in this study lend themselves to a closed-form solution, the transition 
matrix of Equation 2-35 can be used for propagation of the 

covariance matrix from one measurement to the next as follows, where a 
new notation is adopted so that subscripts k and k-1 denote the times 
of measurement (t^^ and : 


t 



0, P 

k,k-l k-1 


+ s °k 4 


(2-87) 


where 


- updated covariance matrix at time tj^_j^ 

I 

- covariance matrix at tj^ before updating 

- are defined in Equations 2-32 and 2-lS, 


The updating equation at time tj^ is given by: 


P 


k 


^ - ■'k «fc '■k 


( 2 - 88 ) 


where 


Kj^ is the Kalman gain matrix defined by 

■‘k = 4 4 t'’k + "k '’k <2-831 

Recursive solution of Equations 2-87, 2-88, and 2-89 provides the 
filtered value of the covariance matrix at any desired time. 



2.7 CONSIDER STATE EQUATIONS 


2.7-1 Introduction 


Optimal performance of the Kalman filter depends on accurate 
description of the statistics of the measurements involved. This im- 
plies that the estimator should accurately model all real world states. 
It may develop that some states can be neglected because they are 
truly negligible. In other cases the statistics of the states may be 
poorly known or the states may incorporate non-zero biases. Here the 
attempt to incorporate poorly modeled states in the estimator may result 
in poor performance, since the estimator is trying to fit the data to 
an incorrect model. 


A useful approach is to have the estimator "consider" the 
effect of poorly modeled states (but not to make estimates of them) 
when estimating the "desired" states. An important advantage of this 
approach is that the number of estimated states is not expanded while 
the effects of "consider" states are still taken into account. 

2.7.2 Consider State Equations 

The equations for the estimation problem with "consider" states 
are as follows: 

Define 



where 

X is the twelve dimensional vector defined by 
Equation 2-12 

is the " consider" portion of the new state vector 
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The covariance matrix for is: 




— _ 


p 

c 

= 

T 

C W 

UP J 

(2-91) 

where p 

= E 


(2-92) 

W 

= E 

(x X '^) 
— c — c 

(2-93) 

C 

P 

= E 

T 

(2-94) 

2 . 7. 3 Consider 

State 

Dynamics 



Since a circular orbit with a nominal attitude history is ass’jaaed 
in this study, a transition matrix can be used here for propagation of 
the covariance matrix. The covariance propagation equation for the 
state ^ is 

K = *T + °T OtG? (2-95) 


twhere the subscripts k and k-1 , to denote times of measurement , are 
omitted for simplicity) and where 


- covariance matrix at time after update, 

{before propagation) 


- covariance matrix at time tj^ after propagation, 
(before update) 


Q 


T 


Q 

0 


G 


T 


G 

0 



0 

G 


(2-96) 


(2-97) 
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Q 


c 


white noise matrix associated with the consider 
states 


- identity matrix (12 x 12 here) 


0 


(2-98 ) 


where <t> is the consider state transition matrix, 
c 

Substituting Equations 2-91 and 2-97 into 2-95, we obtain for the 
propagated covariance matrix the following : 


where 



' T T 

P = $ P ({> + gqg 
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0 c. 
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T 

cf, w 

c c 


T 

t G Q G 
c c c 
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( 2 - 100 ) 

( 2 - 101 ) 


( 2 - 102 ) 


For the present study the consider states are taken to be the 
ephemeris uncertainties in position and velocity associated with each 
of the two TDRSS satellites stationed over the equator. 

The consider state transition matrix, 4*^, is then given by: 


where and are the six-by-six transition matrices associated 

with the ephemeris uncertainties for each of the two TDRSS satellites. 
These matrices are identical in form to the transition matrix given in 
Equation 2-34 for the EOS ephemeris state. In this case the term of 
Equation 2-34 is the orbital rate of TDRSS. 


cl 


0 

^c2 


(2-103) 
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2.7.4 Consider State Measurement Equations 


The measurement equation incorporating the consider states is : 


or 


z 


[H j H^] 



At n 


Z = Ht Xt + n 

where 

Hj = [H I H^] 


(2-104) 


(2-105) 


(2-106) 


and 


H is the sensitivity of the measurement to the state, x. 

is the sensitivity of the measurement to the consider 

state, X . 

-c 

2 . 7. 5 Consider State Filter Equations 

The Kalman filter equations for the covariance are as fallows 
where the updating equation at time, tj^ , is given by: 


P 

c 



(2-107) 


where P 


c 


and P^ are the covariance matrices before and after updating. 


and 


T 

^ «T 


[Ht Pe 


At ^ ta At 


(2-108) 

( 2 - 100 ) 
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Here there is no interest in updating the consider state, x , but 

~c 

only the state, x. Therefore, only those rows of which pertain to 
X are calculated. This gain is called K . That is, choose K 
such that: „ _ 


K 


T 


K 

X 

0 


( 2 - 110 ) 


Combining Equations 2-?l and 2-106 yields: 


‘"t ♦ % «c''> K h/ %■’ 


( 2 - 111 ) 


where 

S ^c ^^t'^ ^ Hp' + HCp t H^WH^ + Rp (2-112) 


The updating equations for P, and W can be shown to be 

P = P’ “ (Eip' + H^Cp^) , 

= C ' - K (Hc' + H w' } , 

P P X p c 


W = W 


(2-113) 

(2-114) 

(2-115) 


It should be noted that the update in Equation 2-114 should be 
carried out even if the consider state sensitivity is zero, such 
as would be the case when landmark measurements are used. 


2.8 ERROR STUDY RESULTS 


2.8.1 Introduction 

The results presented in this section give performance data on 

the landmark, range and range-rate update configuration. The previous 

( 1 ) 

report, Interim Technical Report No. 2 , had presented results on 

the landmark-star configuration. 

The present study assumes radar range and range-rate measurements 
made by either (1) , two ground tracking stations (located at Merritt 
Island, Florida and Goldstone, California), or (2), two TDRSS^^^ geosyn- 
chronous navigational satellites located over the equator. Most of the 
results are generated for landmark observation case 2, which includes two 
passes over the continental USA. (See Figure 2-4) . The material presented 
in Section 2.8.^ shows that additional passes over the USA and Alaska have 
a relatively small effect on reducing estimation uncertainties. Figure 2-5 
illustrates the four-pass case and Table 2-1 briefly describes all of the 
cases used in the study. 

The landmark location data is generated artifically, since no 
realistic data was conveniently available nor was considered necessary 
for this type of study. Most performance results are generated using 
only two landmark updates per pass over the continental USA. The land- 
mark location data is obtained at points within the USA by selecting 
random values of the multi-spectral scanner (MSS) beam angle- The 
maximum scan beam angle is + 4.8° corresponding to + 45 miles on the 
ground. (See Figure 2-1) . The landmarks are located near the north 
and south USA borders for the given pass. 

For the ground tracking stati''^ s -f- _ . h ri ons , five simultaneous 
range and range-rate updates per puss are assumed. The first pass is 
tracked by Merritt Island and the second by Goldstone. Minimum allow- 
able radar beam elevation is 5 degrees. For the simulations with TDRSS 
tracking, eighteen simultaneous range and range-rate measurements by 
both TDRSS stations are assumed. The TDRSS satellites are located at 
41 and 171 degree.s west longitude over the equator. More TDRSS measure- 
ments are possible because of the greater coverage of EOS orbits from 
the vantage point of the TDRSS satellites. 
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TABLE 2-1 

LANDMARK OBSERVATION CASES 


Case 

Number 

Pass 

North Latitude (deg) | 

Reaion* 

Number 

At: Start 

At End 

1 

1 

50 

30 

USA 

2 

1 

50 

30 

USA 


2 

50 

30 

1 

3 

1 

50 

45 

USA 


2 

50 

30 

1 


3 

50 

45 

Y 

4 

1 

50 

30 

USA 


2 

50 

30 

1 


3 

65 

60 

Alaska 


4 

65 

60 

\ 

5 

1 

50 

45 

USA 


2 

50 

30 

1 


3 

50 

45 

y 


4 

65 

60 

Alaska 


5 

65 

60 



*USA denotes continental USA. 


2-36 





































For the dotnparative simulations using a star tracker, only one 

6 

type of star tracker is assigned. This is th" body fixed 8 square 
FOV tracker which electronically tracks a star as it passes through 
the FOV. For the relevant runs the stars are chosen at equal orbital 
intervals ( 20 stars per orbit) with randomized locations within the 
tracker’s field of view. 

The performance results presented here are generally the uncer- 
tainties in the Kalman filtered estimates of spacecraft attitude, gyro 
bias drift and ocbital ephemeris after processing landmark, range and 
range-rate measurements for two passes. The results are the square 
roots of the principal diagonal terms of the 12 x 12 covariance matrix. 
They provide a good Ptatistical indication of performance in estimating 
the state parameters. 

2.8*2 Nominal Values Used in Study 

Unless otherwise stated, the nominal values of the error sources 
Slid parameters used to generate the performance results are as given 
here. Although some of the values may not represent the best or latest 
estimates, it is felt that they are satisfactory for the present study. 

The landmark measurement uncertainty used in this study (15 
meters) is a composite number representing both uncertainty in the 
landmark location in Earth coorc^inates and uncertainty in the multi- 
spectral scanner measurement. Two assumptions are made concerning the 
landmark. First it is assumed to have no associated recognition pro^ 
cess errors, and secondly it is assumed to be essentially a point 
source. These requirements are met by the searchlight landmark des- 
cribed in Reference 1. There are four basic uncertainty sources in 
this measurement, namely: 

1) artificial landmark surveying errors, 

2) earth pole wander and rotation rate variations, 

3) MSS resolution, 

4) scan rate and timing errors. 
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The uncertainty numbers associated with each of these sources are some- 
what arbitrary, for example surveying accuracy can be increased in 
steps by applying more and more sophisticated techniques. 


Estimates for the uncertainty in each category are as follows: 


surveying 
earth wobble 
MSS resolution 
scan errors 


10 meters, 

3 meters, 
10 meters, 

4 meters. 


The root sum-of -square value for these uncertainties is 15 meters. 
Further discussion of the sources of these numbers is gix'en in section 
3.3. 


The nominal values for uncertainties in the range and range-rate 
measurements (10 meters and 0.012 meter/ sec) include both biases and 
white noises in the measurements. Because of the limitations imposed 
upon the present study the radar biases are not modeled separately as 
estimatable states or as consider biases, since more extensive analysis 
and simulations xtfould have been required. Although the data given in 
the literature for these noises varies somewhat, it is felt that the 
above values are reasonable selections. 

The nominal value for ground tracking station location uncertainty 
is assumed to be 5 meters along each axis. This figure has been chosen 
because Merritt Island and Goldstone are particularly well surveyed, 
although the nominal location uncertainty for all USES (United-S Band 
System) radars is 10 meters. 

TDRSS Ephemeris uncertainties 

The choice of a suitable TDRSS ephemeris uncertainty covariance 
is a difficult probleim. The only available TDRSS uncertainty covar- 
iance was printed in. a 1971 report on the TDRSS system. This co- 
variance was reprinted in a 1973 NASA report on "Navigation Systems 
Characteristics, Rev However this covariance may have repro- 

duction errors since one of the cross-correlation coefficients is 
greater than unity. Nevertheless this covariance is used in the TDRSS 
simulations after appropriate modification. 


ftEPRODUCBILITY OB' mB 
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Two other TDRSS ephemeris uncertainty covariances are also used 
in the simulations for comparative purposes. Data for the first covar- 
iance was obtained from G5FC and represents a somewhat pessimistic 
estimate. These assume position uncertainties about four times greater 
than for the first covariance. Data for the second covariance antici- 
pates future state-of-the-art improvements in determination of TDRSS 
ephemeris. Here position uncertainties are about 50 times smaller than 
for the first set of uncertainties. 


Gravitational Harmonic Phcertainties 

All covariance runs made for this study include the effects of 
uncertainties in the gravitational harmonic c & S coefficients as well 
as in the gravitational constant* y. A separate study was first under- 
taken where the term-by-term differences between the most recent Smith- 
sonian gravitational models (SEII (1969) and SEIIi (1973)) were treated 
as if equivalent to one sigma uncertainties, and where the gravitational 
constant uncertainty was 0.2 PPM. The effect on the computed EOS orbit 
beyond radar coverage (for both GTS and TDRSS) was determined for each 
term. It was found incidentally that the C and S .coefficients having 
the most significant eftect on EOS orbit were those associated with 
the following n , m pair s; 2,2; 3,3; 5,5; 6,2; 6,5; 6,6; 8,7. By sum- 
ming the ephemeris error results of numerous Monte Carlo runs for which 
the C and S coefficients were randomized, an uncertainty covariance 
matrix was obtained that represented the ephemeris uncertainties at the 
point that radar coverage commences prior to passage over the USA. This 
covariance was then scalerl to match data obtained from GSFC for 5 orbits. 

With ground tracking coverage, the square roots of the uncertainty 
covariance diagonals were as follows: 

(altitude, down range, crosstrack) i 

(2.6, 14,3, 3.3, ,013, ,003, .002) meters, meters/sec 

With TDRSS tracking, because of better coverage, the corresponding 
uncertainties were: 


(3.0, 6.8, 7.8, .007, 003, .003) meters, meters/sec 


For the 12-state covariance simulation studies this ephemeris 
uncertainty covariance matrix was summed with the propagated 12-state 
covariance at the point for each pass when radar coverage commences 
prior to USA passage. 

The nominal values of the error sources and parameters are 
as follows: 


Initial State Uncertainties (lo) 


Attitude (Pitch, Roll, Yaw) 

- 

60 

arcsec (each) 

Gyro Bias Drift 

- 

0.03 

deg/hr (each) 

Ephemeris Position 

Attitude 

- 

20 

meters 

Downrange 

- 

50 

meters 

Cross track 

- 

20 

meters 

Ephemeris Velocity 

Altitude 

- 

0.05 

meter/sec 

Downrange 

- 

0.02 

meter/ sec 

Cross track 

- 

0.02 

raeter/sec 

^ro Error Sources (la) 

Random Drift 

- 

0.01 

deg/hr (white 

Quantization 

- 

0.1 

arcsec 


Landma rk Measurements (la) 


Landmark Position 
Downrange 
Crosstrack 


15 meters 
15 meters 


flo fqv star Tracker 

Measurement Error Cla) 
Field-of-View 
Pointing Direction 
Star Distribution 


5 arcsec/axis 

8x8 degrees. 

towards zenith 

star randomly selected 

in FOV after each 18 

degrees of orbital motion 
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Range and Range Rate Measurements (la) 


Slant Range 10 meters 

Slant Range Rate 0.012 meter/sec 

GTS Location 5 meters (each) 

TDRSS Ephemeris (Position & Velocity; 

Altitude, Downrange, Crosstrack) 

1971 Report 

(15, 276, 128; .020, .001, .009) meter ,raeter/sec 
GSFC 

(150, 1000, 500; .073, .011, .036) meter, meter/sec 
Future State-of-the-Art 

{5,5, 5; . 00037, .00037, .00037) meter , meter/sec 

Gravitational Coefficient Uncertainties (lo) 

Gravitational Constant, y 0.2 PPM 

C&S Coefficients (See previous comments) 

Normal Conditions for Observation Runs 

Landmark Observation Case 2 (Includes 2 passes over 
continental USA) 

2 Landmark updates per pass 

5 Range and Range Rate updates (simultaneous) by 
Merritt Island on first pass and by Goldstone on 
second pass 

Run starts at equator at 98° east longitude on ascending 
node and passes over north polar region before making 
first pass over USA 

Run ends just south of last landmark update on second 
pass over USA 

2.8.3 Comparative Results yith Different Measurement Configurations 

Before presenting the re alts obtained with different measure 
ment configurations, it will be useful to show the effects of using 
single type of measurement for updating the estimation uncertainties 
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Table 2- 2 gives performance data after processing with a single type of 
measurement. All data presented includes the effects of two passes over 
the USA. 

For comparison the state uncertainties for the case where no 
measurements are made are presented in the first row of Table 2-2. In 
this ease the covariance matrix of uncertainties is simply propagated 
without update with any measurements. The increase in attitude un- 
certainty is due to the uncertainty in gyro bias drift, while the in- 
crease in satellite altitude and downrange position uncertainty is due 
to the effect of initial altitude uncertainties on the computed gravity 
feedback . 

The next row shows the effect of star sightings with the 8* FOV 
star tracker where it is assumed that 20 stars per orbit could be seen 
with the tracker. The results show a very significant decrease in the 
attitude and gyro bias drift uncertainties. However the ephemeris un- 
certainties were not improved at all. This is because star sightings 
do not provide any information on satellite position. 

The third row of Table 2-2 gives the effects of measurements of 
known landmarks where two landmark updates per pass over the continental 
USA were used. The results show only a moderate decrease in attitude 
and drift uncertainties and a relatively slight decrease in ephemeris 
uncertainties. Note that the pitch attitude and ephemeris uncertainties 
are greater than their initial values. 

The next row shows the effect of range and range rate measure- 
ments by ground bracking stations (Merritt Island and Goldstone) where 
it was assumed that the range and the tange rate measurements were made 
simultaneously. The results contrast sharply with those for star sight- 
ings. Here there is a very significant decrease in ephemeris uncertain- 
ties, while there is no improvement in attitude and gyro bias drift un- 
certainties . The last row also shows the effect of range and range rate 
measurements taken from the planned geosynchronous Tracking and Data 
Relay Satellite System (TDRSS) . Although more measurements per pass 
(18) were made (because of better converage of the EOS orbit by TDRSS) , 
the performance is markedly inferior to that with ground trackii.g. This 
is primarily due to the relatively larce initial ephemeris uncertain- 
ties assumed for the TDRSS satellites. This will be discussed in the 
next subsection. 
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TABLE 2-2 

KALMAN FILTER PERFORMANCE WITH 
SINGLE TYPES OF MEASUREMENTS 



: 



State 

Estimation Uncertainties (la) 



Ml 

Type 

of 

Attitude 

(arcsec) 

Gyro Bias 
(10”3 de£ 

Drift 

j/hr) 

Position 
(meters ) 

Velocity 

(meters/sec) 

Measurement 

“pitch] 

Roll 

Yaw 

X 

Y 


Alt. 

Range 


Alt. 

Range j Track j 

Initial State Uncertainties 

. 











60 

60 

60 

30 

30 

30 

20 

50 

20 

.05 

.02 

.02 

No 

Measurements 

273 

83 

83 

30 

30 

30 

101 

745 

20.5 

.749 

.075 

.020 

Stars 

2.0 

1.8 

14.9 

15 

0.4 

1.0 

101 

745 

20.5 

.749 

.075 

.020 

Landmarks 

94 

4.7 

26.8 

24 

12 

22 

57 

460 

20.0 

.482 

.039 

.020 

Range 

and Range Rate 
(GTS) 

273 

83 

83 

30 

30 

30 

2.8 

4.6 

7.8 

.012 

.003 

.020 

Range 

and Range Rate 
(TDRSS) 

273 

83 

83 

30 

30 

30 

10.1 

44.4 

16.8 

.040 

.010 

.016 


Notes: Nominal numlser of measurements per pass or per orbit 

are assumed for each measurement type. See Section 2.7.2 
GTS - Ground Tracking Station 



























Table 2-3 shows the performance for various measurement configura- 
tions. The first row shows the results when using both star sightings 
and known landmark measurements. The decrease in attitude uncertainty 
is due almost entirely to star information, since it is almost the 
same as with stars alone. However the decrease in attitude uncertainty 
enables the landmark information to be used more effectively to reduce 
the ephemeris uncertainty. The next three rows give the state uncer- 
tainties that result from a measurement configuration that combines 
known landmarks and range and/or rarge rate ground tracking. The de- 
crease in ephemeris uncertainty is essentially the same with range and 
range-rate updates taken alone, as compared with the configuration where 
landmarks are included in the measurements. However the existence of 
range and range rate information enables the landmark measurements to 
be very effectively used to reduce the pitch and roll attitude uncer- 
tainties as well as the Y bias drift uncertainty. The other two rows 
show that range-rate measurements are more effective than range measure- 
ments in decreasing ephemeris uncertainties. 

The next to last row of Table 2-3 gives the results of using stars 
and range and range-rate information with no landmarks. Here the de- 
crease in attitude uncertainty is due entirely to star sightings, while 
the decrease in ephemeris uncertainty is due entirely to range and 
range rate information. The addition of landmark measurements to this 
configuration decreases attitude uncertainties slightly to the lowest 
uncertainty level of any combination considered. 

2.8.4 Comparative Results with Ground Tracking and TDRSS Satellites 

Table 2-4 compares the performance of ground and TDRSS tracking 
for the cases with and without uncertainties in the positions of the 
ground stations or TDRSS. Landmark measurements were also used in all 
of these cases. The ground tracking stations had a nominal location un- 
certainty of 5 meters per axis. The TDRSS satellites were assumed to 
have initial ephemeris uncertainties of 15, 276, 128;. 020, .001, .009 
meters and meters/sec in altitude, downrange and crosstrack. The 
associated covariance matrix is from a 1971 report. (See Section 
2.8.2), For Table 2-4 the TDRSS ephemeris uncertainties are modeled as 
consider states. 


TABLE 2-3 

KALMAN FILTER PERFORMANCE WITH 
DIFFERENT TYPES OF MEASUREMENTS 


Types 

of 

Measurements 


State Estimation Uncertainties (la) 


Attituae 

{arcsec) 


Pitchl Roll I Yaw j X 


Initial State Uncertainties 


60 6 


Gyro Bias Drift 
(10“3 deg/hr) 


Z 


Position 

(meters) 


30 30 30 


20 50 


Landmarks 

and 

Stars 


Landmarks 
& Range 


Landmarks & 
Range Rate 


Landmarks & 
Range & 
Range Rate 


Stars & Range 
& Range Rate 


Landmarks & 
Stars & Range 
& Range Rate 



.6 22 2.8 



Velocity 

(meters/sec) 

1 Alt. 

Range Track 

.05 

.02 

.02 

.030 

.015 

.02 

.016 

.007 

.020 

.014 

.003 

.020 

.012 

.003 

.020 



Notes: Nominal number of measurements per pass or per orbit are assumed for each measurement 

type. See Section 2.7.2. Ground tracking stations used for range and range rate 
measurements . 


a * li I I'l' r r Villi' it *iiTi AjWsrlAl&hl 























































TABLE 2-4 

KALMAN FILTER PERFORMANCE WITH 
GROUND AND TDRSS TRACKING STATION MEASUREMENTS 


Tracking Station 

State Estimation Uncertainties (lo) 

Ephemeris 

Attitude (arcsec) 

Position (meters) 

Uncertainties 

Pitch 

Roll 

Yaw 

Alt. 

Range 

Track 

Initial Uncertainties Are: 

60 

60 

60 

20 

50 

20 

Ground Tracking Stations 







Zero Locat.Uncert. 

2.4 

2.5 

26.8 

2.1 

2.8 

6.6 

Nom . Local . Oncer t . 

2.5 

2.7 

26.8 

2.8 

4.6 

7.8 

TDRSS Satellite Stations 






■jH 

Zero Ephem.Uncert. 

3.2 

2.3 

26.6 

4,8 

11.2 

IBH 

Nom. Ephem. Uncert . 

9.2 

4.1 

26.7 

9.7 

44.3 

16.8 


Note: All data obtained with nominal number of range 

and range rate and landmark measurements. 

TDRSS ephemeris uncertainties modelled as consider states. 











Table 2-4 shows that with zero ground station and TDRSS location 
uncertainties, the position estimation uncertainties with TDRSS meas- 
urements compare favorably with ground tracking results, particularly 
in crosstrack. However with nominal ground station and TDRSS location 
uncertainties the position estimation uncertainties are markedly worse 
with TDRSS tracking. Note for example that the ^^vnrange position 
uncertainty is only slightly less than the initial uncertainty. 

TDRSS Ephemeris Uncertainties Modeled as Consider States 

The TDRSS ephemeris uncertainties are associated with a state 
that is highly time-correlated over the simulation time intervals of 
this study. Since the errors in this state will not undergo very much 
change during those time intervals, it is appropriate that they be 
modeled as consider state biases rather than as zero bias white 
noises. With this type of modeling the Kalman filter estimator 
‘‘considers" the effect of imperfectly modeled states (but does not 
measure or estimate them) on the estimation of the "desired" states. 
Since the TDRSS ephemeris uncertainties may have non-zero means rel- 
ative to the data intervals studied, the treatment of these uncer- 
tainties as "consider" states implies better modeling of the real 
world data. The analysis and equations for "consider" state pro- 
pagation and update are presented in Section 2-7. 

While it is better to model the TDRSS ephemeris uncertainties as 
consider stater^ rather than as white noises, it is nevertheless of 
some interest to investigate Kalman filter performance with both mod- 
eling types. Table 2-5 compares estimation performance for the two 
types of modeling with three sets of initial TDRSS ephemeris uncer- 
tainties. These were given in Section 2.8.2 and are repeated here 
for the readers convenience. The three sets of initial ephemeris un- 
certainties are: (1) an adverse estimate obtained from 6SFC, (2) a 
conservative estimate first given in a 1971 report, and (3) a future 
state-of-the-art estimate now being used in other unrelated CSDL 
studies. The one sigma ephemeris uncertainties in meters and meters 
per second for the three sets are: 
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TABLE 2-5 

KALMAN FILTER PERPORMANCE WITH 
DIFFERENT TDRSS EPHEMERIS ONGERTAINTIES 


Tracking Station 
Ephemeris 
Uncertainties Clo’) 


State 

Estimation 

Uncertainties (la) 

~ 1 

Attitude ( ar csee ) 

Position (meters) | 

Pitch 

Roll 

Yaw 

Alt- J 

Ra. . je 

Track 

Initial Uncertainties Are: 








€0 

60 

60 

20 

50 

20 

0 

3.2 

2.3 

26.6 

4.8 

11.2 

3.1 

With GSFC Ephemeris Uncertainties Modelled as: 





Noises 

23.7 

4.2 

27.0 

18.9 

114.4 

17.7 

Consider States 

39.6 

4.6 

27.1 

27.2 

195.1 

19-9 

[with *71 Shuttle Study Ephemeris Uncertainties 

Modelled as : 



Noises 

6.5 

2.8 

26.7 

9.7 

30.0 

8.4 

Consider States 

9.2 

4.1 

26.7 

9.7 

44.3 

16.8 

With Future State-of-the-Art Ephemeris Uncertainties Modelled as: 



Noises 

3.3 

2.3 

26.6 

4.9 

11-9 

3-3 

Consider States 

3.6 

2.3 

26.6 

5.1 

13-7 

3.7 


Note: All data obtained with nominal number of range 

and range rate and landmark measurements. 










Alt. 

Dnrng . 

Crtk. 

^alt 

'^dnrng 

'^crtk 

GSFC conserv. 

1) 

150, 

1000, 

500; 

.073, 

.011, 

.036 

1971 rpt 

2) 

15, 

276, 

128; 

.020, 

.0011, 

.0025 

Future est. 

3) 

5, 

5, 

5; 

.00037, 

.00037, 

.00037 


Table 2-5 shows that the Kalman filter performance is consistent- 
ly worse when the TDRSS ephemeris uncertainties are modelled as "con- 
sider” states than when modelled as white noise. This is to be antici- 
pated, since use of "consider” states represents a more realistic model- 
ling of possible biases. Table 2-5 also shows that with consider state 
modelling, the Kalman filter performance for ephemeris uncertainty esti- 
mation is only satisfactory with future state-of-the-art TDRSS ephemeris 
uncertainties. With 1971 Shuttle study uncertainties the downrange 
estimation uncertainty is only slightly less than the initial uncertain- 
ty. 


2,8,5 Filter Performance for the Five Landmark Observation Cases 


Table 2-6 presents the Kalman filter estimation uncertainties 
at the end of the last pass for the five different landmark obser- 
vation cases* (See Table 2-.1) . Numerically these cases correspond to 
the number of passes over the continental USA, and if Case 4 or 5, 
over Alaska aa well. To track the last two passes a third USES radar 
tracking station located at Fairbanks, Alaska, is used for range and 
range rate measurements. 

Table 2-6 shows that most of the improvement in performance is 
accomplished by the end of the second pass. Only slight improvement 
results as the number of passes is increased from two to five. In 
going from Case 2 to 3 it is noted that the epheraerls uncertainties 
did not decrease as one might expect. For example, the downrange un- 
certainty is 6.6 meters for Case 3, while it is 4.6 meters for Case 2. 
The reason for this irregularity lies in the choice of EOS orbital 
paths for the different cases. To make possible three observational 
passes over the continental USA (Case 3) , a somewhat different inertial 
orbital path from that for Case 2 had to be used. Because of the 
different relationships of the EOS orbital paths with respect to the 
ground tracking stations for Cases 3 and 5, as opposed to those for 
Cases 2 and 4, the estimation uncertainties are moderately different 
from what they would have been if the same inertial orbital path had 
been used for all cases. 

2.8.6 Kalman Filter Performance as a Function of Orbit Angle 
Position Estimation Performance 

Figures 2-6 and 2-7 show graphically the effect of propagation 
and update of the position estimation uncertainties for Cases 2 and 4 
when using landmark measurements and ground tracking range and range 
rate measurements. For comparison, Figure 2-8 shows the position un- 
certainties for Case 2 when using star sightings and the same ground 
tracking mesurements. 
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TABLE 2-6 


KALMAN FILTER PERFORMANCE FOR 
DIFFERENT LANDMARK OBSERVATION CASES 


Landmark 

State Estimation lancer taint ies 

(la) 

Observation 

Attitude (arcsec) 

Position (meters) 

Case 

Pitch 

Roll 

Yaw 

Alt, 

Range 

Track 

Initial Uncertainties Are; 








60 

60 

60 

20 

50 

20 

1 

IQI 

5.2 

34.3 

12.3 

9.7 

19.2 

2 (nom) 

2.5 

2.7 

26.8 

2.8 

4.6 

7,8 

3 


1.9 

21.2 

3.2 

6.6 

wm 

4 

B 

1.9 

18.4 

2.3 

3.6 

6.2 

5 

1.9 

1.6 

16.9 

B 

B 

4.8 


Notes: 2 Landmark Updates/Pass 

5 Range and Range Rate Updates/Pass by Ground Tracking 
































M 

i 

LH 

(NJ 


P SO 

§ 2 

£> O 

Sp 

a 


" H 

aj 



Figure "^-G Kalman Filter Position Estimation Performance 
for Lanclnark Observation Case 2 
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POSITION UNCERTAINTY (m) 


NOTES : 

2 Landmark Updates/Pass 

5 Range 6 Range Rate Opdates/Pass 

Orbit Angle Measured from Ascending Node 


■Downrange 


•Over USA 


•Altitude 


Cross* 

track- 




Over 

Alaska 


n 
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Over 

Alaska 


ORBIT ANGLE (deg) " 

Figure 2-7 Kalman Filter Position Estimation Performance 
for Landmark Observation Case 4 
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Figures 2-6 and 2-7 show how effectively range and range rate 
measurements reduce the altitude and downrange position uncertainties 
on the first update of each pass for Cases 2 and 4. Subsequent updates 
on i-he Scune pass only reduce the uncertainties slightly. The marked 
increase in downrange uncertainty #and to a lesser extent the altitude 
position uncertainty, after the first USA pass is due to the presence 
of relatively large velocity estimation uncertainties. The operation 
of the Kalman filter on the velocity uncertainties is shown in Figure 
2-7 by the marked reduction in peak downrange and altitude position 
uncertainties resulting after the second and third passes. The cross- 
track uncertainty is more cyclical than the other uncertainties. This 
uncertainty is progressively reduced after each pass. 

Figure 2-8 shows graphically for Case 2 the effect of using star 
sightings with the same range and range rate measurements. It is seen 
here that the reduction in the downrange and altitude uncertainties at 
the start of each pass over the USA is not as dramatic as was the case 
in Figures 2-6 and 2-7 where landmarks were used. 

Attitude Es t imat ion Per f ormanc e 


Figures 2-9 and 2-10 show graphically the attitude estimation 
performance when using landmark, range, and range rate measurements. 
Figure 2-11 does the same with landmark and star measurements. 

Note in Figures 2-9 and 2-10 that significant reductions occur 
in the attitude uncertainties (particularly those of pitch and roll) 
during the passes over the USA and Alaska. In these cases, the only 
information provided on attitude was that obtained from the landmark 
measurements during the passes over the USA and Alaska. Two landmark 
measurements were made on each pass. Most of the reduction in attitude 
uncertainty during a pass occurs during the first landmark update of 
each pass. In addition, it is seen that most of the reduction occurs 
in pitch and roll since the landmark is always near the local vertical. 
Note that the performance in Case 2 (Figure 2-9) is essentially the 
same as that of Case 4 (Figure 2-10) during the first two passes since 
Case 4 is essentially an extension of Case 2. Also note in Figure 2-10 
that the uncertainty in pitch remains low and does not grow as rapidly 
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Figure 2-10 Kalman Filter Attitude Estimation Performance 
for Landmark Observation Case 4 



after the second pass. This is due to the fact that the uncertainty 
of the Y <pitch) gyro bias drift was greatly improved by the updates 
of the second pass. 

In Figures 2-9 and 2-10 it is seen that the roll and yaw uncer- 
tainties become oscillatory after the first pass. Note that the fre- 
quency in yaw is twice that in roll, and that there is a definite phase 
relationship between the two. The reason for the oscillatory behavior 
is due to a complex relationship between roll, yaw, roll gyro bias 
drift, and yaw gyro bias drift, and the fact that updates are made only 
during the passes over the continental USA and Alaska, which occur at 
somewhat the same portion of the satellite orbit. 

In dramatic contrast to Figure 2-9 and 2-10, Figure 2-11, with 
star sightings in combination with landmark measurements, shows a very 
different performance in attitude estimation. This is primarily due 
to the fact that many more updates are made in attitude (20 star sight- 
ings per orbit plus 2 lancmark measurements per pass) . Note that the 
frequent updates of attitude throughout the orbit causes the uncertain- 
ties of all attitude components to reach almost steady state values 
after the first orbit. 

2*8.7 Performance Sensitivity to System Errors and Parameters 

In this section the sensitivity of Kalman filter performance to 
various system uncertainties and parameters is inve&f:igated. These 

■k 

include the following: 

1) Initial ephemeris uncertainty 

2) Initial gyro bias drift uncertainty 

3) Ground tracking station location uncertainty 

4) Range measurement noise 

5) Range rate measurement noise 

6) Landmark position uncertainty 

7) Number of range updates per pass 

8) Number of range-rate updates per pass 

9) Number of range and range-rate updates per pass 

♦No data is presented on the sensitivity to initial spacecraft 
attitude uncertainties since the performance in estimating the 
ephemeris and attitude was virtually unaffected by values of these 
uncertainties up to 15 times nominal. 



Figure 2-11 Kalman Filter Attitude Estimation Performance 

for Landmark Observation Case 2 with Star Sightings 




10) Number of landmark updates per pass 

11) Gyro random drift and quantization 

12) Gravity harmonic uncertainties 

These are considered separately in the following paragraphs . 

Sensiti^rity to Initial Eplhiemeris Uncertainty 

Figure 2-12 shows the effect of variation of the initial ephemeris 
uncertainty (all six components varied simultaneously) on position and 
attitude estimation uncertainties. This figure demonstrates that even 
when the initial ephemeris uncertainty is ten times its nominal value# 
the altitude and crosstrack position uncertainties are only slightly 
affected, while the downrange uncertainty is moderately increased. As 
might be expected the effect on the attitude uncertainties is neglig- 
ible. 


Figure 2-13 shows that only the yaw attitude uncertainty is 
affected by increases in the initial gyro bias drift uncertainty about 
all axes. In this case it is seen that the position estimation un- 
certainties remain unaffected. 

Sensitivity to Ground Tracking Station Location Uncertainty 

Figure 2-14 shows that the performance in estimating spacecraft 
position is strongly affected by an increase in the tracking station 
location uncertainty (assumed equal in all dimensions) . For example, 
an increase in the location uncextainty from 5 (nominal) to 10 meters 
increases the downrange estimation uncertainty from 4.6 to 6.9 meters. 
This is to be expected since the location uncertainty effectively adds 
to the range measurement error. 

Sensitivity to Range Measurement Noise 

Figure 2-15 shows the effect of variation in the range measure- 
ment noise while maintaining nominal values for the errors in range 
rate and landmark measurements. It is seen that as the range measure- 
ment noise is increased# the performance in estimating spacecraft 
position levels off at values which happen to be those primarily obtained 
with range rate measurements alone. Mote that there is little to be 
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Figure 2-12 Sensitivity to Initial Ephemeris Oncertainty 









Figure 2-13 Sensitivity to Initial Gyro Bias Drift Uncertainty 
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Figure 2-15 Sensitivity to Range Measurement Noise 






ir^-.i-r ■ . , I, ..Iii^. ...^ I I II P H I I f ■ WFW in IH r.... . i 


gained from range measurements when the one-sigma noise in these 
measurements exceeds about twice the nominal value of IQ meters. 

Sensitivity to Range Rate Measurement Noise 

Figure 2-16 shows the effect of variation in the range rate 
measurement noise while maintaining all other conditions at nominal. 

In contrast to the results just shown for range measurement noise in 
Figure 2-15, it is seen that the sensitivity to range rate noise ex- 
tends over a much larger range of variation in this parameter. This 
is also seen to be the case in Table 2-3. 

Sensitivity to Landmark Position uncertainty 

Figure 2-17 shows the sensitivity of performance to landmark 
position uncertainty when landmark measurements are used in combination 
with ground tracking measurements. It is seen that the performance in 
estimating attitude is strongly affected by a variation in the land- 
mark position uncertainty, while the performance in estimating space- 
craft position is essentially unaffected. This is as expected since 
the landmark measurements provide the only data on attitude, while the 
ground tracking measurements, with the assumed high accuracies, dominate 
the performance in estimating spacecraft position. If landmark measure- 
ments had been used with star sightings instead of ground tracking 
measurements, the reverse situation would have occurred, where the star 
sightings dominate the attitude performance and the landmark measure- 
ments provide the only information on spacecraft positiion. (See 
Figure 2-18 of Reference 8) . In either ease, it is important to note 
from Section 2.7.2 that the landmark position uncertainty is treated 
in this covariance study as being the overall error associated with a 
landmark measurement. 

Sensitivity to Number of Range Updates per Pass 

Figure 2-18 shows that approximately 15 range updates per pass 
are required to reach asymptotic levels in position estimation per- 
formance when no range rate measurements are included. Although the 
conditions under which this data was generated are not considered to 
be quite nominal since the nominal range rate measurements were not 
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included, the data at least shows the effect of varying the number of 
range measurements in combination with the nominal number of landmark 
measurements. Also note in Figure 2-18 that only two or three range 
updates are required per pass in order to achieve minimum attitude 
estimation uncertainties. 

Sensitivity to Number of Range Rate Updates per Pass 

Figure 2-19 indicates that minimum position estimation uncer- 
tainties are achieved after about 5 range rate updates per pass when 
no range measurements are included. Note that it takes fewer range 
rate updates per pass, as compared to range updates per pass, to 
achieve minimum estimation uncertainties. Also note that only two or 
three range rate updates per pass are required to reduce the attitude 
uncertainties to minimum levels. 

Sensitivity to Number of Range and Range Rate Updates per Pass 

Shown in Figure 2-20 is the sensitivity of performance to the 
number of updates per pass with both range and range rate measurements. 
Here it is seen that the affect on performance is not much different 
from that in Figure 2-19 where no range measurements were made. In 
addition, it is seen that five range and range rate measurements (the 
nominal number used in this study) yield results which are fairly close 
to what would be obtained if more range and range rate measurements 
were used during a pass. 

Sensitivity to Number of Landmark Updates per Pass 

As was true in Figure 2-17 for sensitivity to landmark position 
uncertainty. Figure 2-21 shows that the attitude estimation uncer- 
tainties are reduced by increasing the number of landmark updates per 
pass, while the uncertainties in estimating spacecraft position remain 
unaffected. Note that only two landmark updates per pass are required 
to bring the attitude uncertainties close to their minimum levels. 

Again it is important to note that the results in Figure 2-21 (like 
those in Figure 2-17) were generated for the case of no star sightings. 
If landmark measurements had been used with star sightings instead of 
range and range rate measurements, the performance in estimating 
spacecraft position would have been affected to some extent by the 
number of landmark updates per pass. (See Figure 2-13 of Reference 8). 
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Figure 2-20 Sensitivity to Number of Range & 
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Sepsitivity to Gyro Random Drift and Quantization 


Performanee results were generated with various values of gyro 
random drift up to six times the nominal value of 0.01 degree per hour. 
For this amount of variation the results showed that random drift had 
no affect on the position estimation performance and only a slight in- 
fluence on attitude estimation performanoe. For example, a six sigma 
random drift on each gyro resulted in pitch, roll and yaw uncertainties 
of 2.6, 3-1, and 27.1 arcsecs, respectively, while nominal values of 
random drift resulted in 2.5, 2.7, and 26. S arcsecs. 

Studies of sensitivity to gyro quantization noise also produced 
similar results. The attitude uncertainties for a six sigma quanti- 
zation noise on each gyro were almost the same as those for a six sigma 
random drift on each gyro. 

Sensitivity to Gravity Harmonic Uncertainties 

In Section 2.1 the gravity harmonic uncertainties were defined as 
being equal to a constant times the difference between the corresponding 
C and S coefficients of the two most recent Smithsonian gravitational 
models (SEIl and SEIII) . This constant (0.33) was chosen so that the 
resulting ephemeris error data would match the results obtained from 
GSFC for five orbits. Inspection of equivalent data from other studies 
indicate that the above definition is probably pessimistic. In Table 
2-7 some performance results are shown for various values of the gravity 
harmonic uncertainties relative to the nominal values defined above. 

This table clearly shows that the indicated changes in the harmonic un- 
certainties have negligible affect on the attitude estimation uncer- 
tainties, while the affect on the position estimation uncertainties is 
moderate . 
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SEMSITIVITY TO GRAVITATIONAL HARMONIC UNCERTAINTIES 
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2. 9 COMCLUSIONS 


In contrast to the previous study where star and landmark 
measurements were used to estimate spacecraft attitude, orbital 
ephemeris, and gyro bits drift, the present study is primarily con- 
cerned with how well the same state parameters can be estimated when 
using landmark measurements in combination range and/or range rate 
measurements with respect to ground tracking stations or the proposed 
Tracking and Data Relay Satellite System (TDRSS) . For purposes of 
comparison, some results of the previous study are included with the 
results of the present effort in Section 2,8. in the previous study 
it was found that landmark measurements alone (See Table 2-2) were not 
sufficient to estimate all twelve state parameters, (At least, this 
was the case for landmarks visible to the raultispectral scanner of 
Landsat or the proposed EOS) . However, by using both star and land- 
mark measurements, it was found that very good performance could be 
obtained in estimating all twelve state parameters (See Table 2-3} . In 
this case, the star sightings were the primary source of data for esti- 
mating spacecraft attitude and gyro bias drift, and the landmark meas- 
urements provided the only information on orbital ephemeris. However, 
when landmark measurements are used only with range and/or range rate 
measurements (See Table 2-3) , the landmark measurements play a reverse 
role in that they provide the only data for estimating spacecraft 
attitude and gyro bias drift, while the range and/or range rate meas- 
urements dictate the performance in estimating orbital ephemeris. In 
comparing the attitude performance for these two cases in Table 2-3, 
it is seen that the performance is noticeably better when using star 
sightings. However, the following differences, which favor the star 
tracker, should be noted: 

1) More star measurements (20 per orbit) were trade than 
landmark measurements (2 per pass) . 

2) The star measurements were uniformly distributed 
throughout each orbit, while the landmark measure- 
ments were restricted to that part of the orbit over 
the USA. 
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Some of the stars were further from the local 
vertical than any of the landmarks due to the 
larger field-of-view of the star tracker (8“ 
square) . Consequently, the i^erformance in 
estimating yaw attitude should be better for 
the star measurements. 

4) The assumed accuracy of the star measurements 
(5 aresecs per axis) is somewhat optimistic for 
the star tracker being considered for this type 
of mission, while the accuracy assumed for the 
landmark measurements may conceivably be more 
realistic. 

However, regardless of the relative merits of star and landmark measure- 
ments in estimating spacecraft attitude, it is seen from the general 
results of this study that landmark measurements are very versatile 
in that they can be used as the primary or backup data source for 
attitude or orbit determination. 

In this study, performance results were generated for landmark 
measurements in combination with range and range rate measurements 
from either TDRSS or ground tracking stations. {See Table 2-4) . In the 
case of TDRSS it was found that the performance was not as good as with 
ground tracking stations primarily because of the larger uncertainties 
adopted for the ephemeris of TDRSS. If these uncertainties were re- 
duced by two orders of magnitude, the resulting performance in esti- 
mating spacecraft attitude and orbital ephemeris would probably be com- 
parable with that obtained with ground tracking. 

From the sensitivity situdies made with landmark measurements, 
combined with range and range rate measurments with respect to ground 
tracking stations, the following conclusions can be made: 

1) Only two observation passes over the USA are required 
for reasonable performance in estimating both attitude 
and orbital ephemeris. 
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It is not necessary to use both range and range 
rate measurements with landmark measurements in 
order to obtain good estimation performance. For 
example, range rate measurements combined with 
landmark measurements do almost as well as when 
range measurements are included. (See Table 2-3) . 

3) Only two landmarks are required on each pass over 
the USA in order to achieve almost steady state 
performance in estimating spacecraft attitude. 

(See Figure 2-21) . 

4) Variation of the landmark position uncertainty and 
the number of landmark updates per pass had no 
affect on orbital ephemeras estimation. Only the 
attitude performance is affected by these para- 
meters. (See Figures 2-17 and 2-21) . 

5) About two updates per pass are required in range or 
range rate in order to achieve near-steady state 
performance in estimating pitch and roll. (See 
Figures 2-18, 2-19, and 2-20), 

6) The performance in yaw is essentially independent 
of range or range rate measurements. 


SECTION 3 


STAR AND iiANDMARK MEASUREMENT EQUATION DEVELOPMENT 
3.1 INTRODUCTION 

In this section the results are presented for the study task 
pertaining to the development of detailed measurement equations for 
the basic onboard navigation measurements. Since the navigation 
system for EOS may rely extensively on computations done at a ground 
processing center, it may be possible to treat the larger state 
vectors associated with a more detailed set of measur'^’nent models. 

Two types of measurements are considered; namely a star direc- 
tion measurement using a star camera device, and a landmark direction 
measurement using a multispectral scanner. The scanner is assumed to 
be one that scans from side to side at a constant rate. 

The previous treatments of the measurements for these two 
sensors assumed that the basic measurements were angles and that 
the associated errors were white noise. In fact the basic measure- 
ments are voltages and time intervals. These are converted to angles 
by means of computations which introduce errors beyond those arising 
in the measurement of the basic quantities. 

If the basic measurements are taken to be angles, as was done 
previously, the errors arising from conversion of basic measurements 
to angles must be accurately represented. h?>'^ever it is simpler and 
more direct to enter the basic measured quantities into the navigation 
estimation equations. This approach may only be tractable when ground 
based computations are considered since use of the raw measurements 
along with more detailed error models implies increased computational 
requirements . 
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3.2 STAR TRACKER MEASUREMENT EQUATION 


3.2.1 Introduction 

Several potential navigation system configurations for EOS 
involve star tracking devices. A leading candidate for this applica- 
tion is a star "camera" device such as the Ball Brothers CT401 star 
tracker. This type of tracker has the attribute that it acquires and 
tracks stars by means of a scanning electron beam, and therefore 
requires no moving mechanical parts. It is a body-mounted device 
whose optical axis is often aligned with the spacecraft negative Zg- 
axis (local zenith) . 

The task addressed here is to derive a measurement model speci- 
fically for a tracker of this type. The goal of the derivation is a 
model that is very general i.e., includes all the major noise types 
and sources. Essential steps in the derivation include formal state- 
ments of the tracker fundamental measurement and the state-vector , 
measurement-vector relationship, derivation of the formal partial 
derivative equations, and definition of the form of the output-measure- 
ment relationship. 

3.2.2 Basic Tracker Operation and Measurement 

The tracker has a square field of view several degrees on a side, 
and the center of the field is pointed at some selected direction in 
spacecraft body coordinates. There are two inodes of operation, namely 
the search and track modes. In the search mode the scanning electron 
beam is driven over a uniform scan pattern by ramp and stepped voltages 
in the directions of the major axes of the field of view. When a star 
signal, represented by a beam current increase, is detected at a se- 
lected threshold value, the associated search voltages are observed. 
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end their values are used to define the star direction . When a star 
has been "found*', it is then tracked (in the track mode) by a re- 
stricted field scan operation, and changes in star position due to 
spacecraft motion are added to the basic location voltages that were 
obtained at acquisition. Thus the basic star tracker measurements are 
a pair of voltages along nominally orthogonal axes . The relationship 
between these voltages and the angles defining the star location in 
the field of view is not, however, a simple proportion because of 
various distortions of the scan geometry. Instead, a functional rela- 
tionship between the basic analog voltages and the star position 
angles is derived by means of calibration against a test star pattern. 
Discussion of this functional relationship is deferred until later. 
Detailed description of the tracker operation can be found in 
Reference 13. 

3.2.3 Formal Measurement Equation 

The relationship between the measurements, state, and noise 
levels in the device is intrinsically nonlinear. It can be written as 

z=z(x,x,v) (3-1) 

— — * — c 

where £ is vector whose elements are the measured quantities, x is 
the spacecraft and tracker state vector, is a vector of unestimated 
variables ^ that affect the state covariance ("consider" para- 
meters) , and V is a vector of white noises. The state x for this 
problem consists mainly of quantities defining the spacecraft attitude 
such as attitude anglc'=^ and gyro bias drift, and estimatable biases in 
the star tracker. 

Assuming that the estimation process to be applied for navigation 
will be linear, we expand Equation 3-1 in a Taylor series keeping only 
first order terms. The formal linearized measurement equation is thus 
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(3-2) 




3z 

3z 

3^ 

3* + 

— C 



In the case of covariance propagation studies, as for example in 
Section 2 of this report, the 3's represent statistical deviations 
of variables from nominal values. In an actual estimation application 
the 3*s represent deviations of the true scate from the estimated 
state. The partial derivative convention used here is the following; 


so that differentiation of column vector z by column vector x yields 
a matrix H (with elements h.-). The matrices H, = 3z/9^, 

A = 3^/9v are here termed measurement "sensitivity" matrices. Thus 
equation 3-2 may be written 

6z “ H6x + H 6x + A6v 

— — G — C ■ 


In the case of the CT401 star tracker the measurement vector ( 2 ) 
is taken to be the scan voltages mentioned previously, namely 


^-(“1 ■ 

\ y/ 

As implied by the nominal boresight axis of the tracker is taken to 
be the z axis. 


There are many potential sources of contributions to the 
measurement error (tSz in Equation 3-2). First, there are possible 
state uncertainties consisting of spacecraft attitude and several 
tracker bias errors that might be estimated. Secondly there are un- 
estimated biases or "consider" biases, ^ and thirdly ., there are 
white noise sources. The major error sources in the star tracker 
measurement are listed below in Table 3-1 along with representative 
one sigma values. 
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TABLE 3-1 


MAJOR ERROR SOURCES, TYPES AND REPRESENTATIVE VALUES* 


Estima table Biases 


Alignment 

- 

5 

arcsec 

Geometric Calibration 

- 

20 

arcsec 

Mechanical Distortion 

- 

5 

arcsec 

Temperature Measurement Bias 

- 

2.5 

arcsec/deg 

Magnetic Field Measurement 

- 

50 

arcsec/gauss 

Consider Biases 

Temperature Fluctuations 

- 

2.5 

arcsec/deg 

Magnetic Field Fluctuations 

- 

50 

arcsec/gauss 

Star Intensity Estimates 

- 

9 

aresec/mag 

White Noise 

Poisson Photon Arrival Rate 

- 

5 

arcsec 

Electronic Noise 

— 

15 

arcsec 


Two bias sources, namely temperature and magnetic field varia- 
tions, can be included in either bias category, depending upon the 


The representative error values are based in part upon 
Reference 13, and in part upon conversations with R,L. Cleavenger 
of Ball Brothers Research Corporation and R. Doxie of MIT's 
Center for Space Research. These errors are very design speci- 
fic and therefore cannot be loosely used for any application. 

Errors in each of the listed categories can be reduced considerably 
by sufficien' application of engineering design effort. 
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details of the tracker data reduction. For example, if temperature 
and axial magnetic field strength are not measured, then deviations 
of these quantities from their calibration values will fluctuate 
with random amplitude on a time scale that is long compared to the 
measurement time. These fluctuations would then be treated as con- 
sider biases. 


3,2.4 Tracker-Inertial Coordinate Relationship 


To evaluate the partial derivatives or seftSitivity mstfioes 
M, A of Equation 3-2 we write a vector equation relating the 

tracker line-of-sight (LOS) to a star, to the inertial directiq;! 
of the star. The ei. *or-free relationship is just 


-B 


(3-4) 


where T^^^ and T^^. are the coordinate rotations from inertial to 
spacecraft body coordinates, as defined in Section 2. S^ is a unit 
vector in inertial coordinates defining the LOS to the star being 
tracked. S_ is the corresponding unit LOS vector to the star in body 
coordinates. To simplify the treatment of tracker errors we define 
another rotational transformation which transforms from tracker 

boresight coordinates to body coordinates. Ip the tracker boresight 
system the star is located in direction S,jv, and Eguation 3-4 is just 


'^Bo'^OI-I "^BT-T 


(3-5) 


The rotational transform T^^ is arbitrary and composed of a pair of 

rotations T^T„ for boresight pointing, and a rotation T for field of 
A B . ^ , 

view orientation. We select the convention of having the star tracker 
2 axis lie along the boresight direction. Thus a star lying along the 
boresight has the following direction in tracker coordinates: 
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This geometry is shown in Pifjure 3-1. The general form for Sjj, is 
obtained by a pair of rotations from tracker coordinates to the star 
direction in tracker coordinates. A rotation about through an 

f ^ 

angle a followed by a rotation about through an angle 3 yields a 
pair of angles in orthogonal planes corresponding to the two electron 
beam angles^ These are related to scan voltages in the tracker. The 
unit vector is thus 

( sin a \ 

- cos oc sin 0 1 . (3“6) 

cos a cos B / 

Errors in the star position measurement are reflected in the 
term on the right in Equation 3-5. 


3.2.5 Measurement Error Relationships 

Errors in alignment of the tracker with the spacecraft axes can 
be expressed as a small angle rotation matrix multiplying 
True values of a and 3 can be represented as the sum of values pre- 
dicted from the voltage-angle calibration relationship, and error 
values. Thus 


= Q + a , 

P e 

3 = . 

P e 

Equation 3-5 can therefore be expressed in terms of the errors as 


Wo A - 


( 3 - 7 ) 


with 



/ sin 

'“p 

+ 






1 - cos 

'“p 

+ 

«e^ 

sin (Bp 

+ 



\ cos 

'“p 

+ 


cos (Sp 

+ 



( 3 - 8 ) 
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If we let a(z) and 6(z) represent the electr^ beam angle func- 
tions of the voltage, then the form of and 3p(z) depends upon 

how one calibrates the tracker distortions. The most straightforward 
calibration is done by fitting a least squares polynomial to simulated 
star test field measurements {Ref .13) . The test field contains a 9x9 
rectangular array of simulated stars, thus the polynomial should not 
be of high order. In Reference 13 first and second order terms are 
used. Following that example we write 


- Aj^(T-T^, B, I-I^) + B, . AjCT-T^, B, r-I^IV^ 

+ A^(T-I^, B, I-I^)v2 . Aj(T-T^, B, . Aj (T-T^, B. i‘I„)V^Vy 


0 = Bj^(T--T^, B, 1 -Iq) + . , . (similar form) 


(3-9) 


The coefficients and are functions of deviations of temperature 

(T) , magnetic field strength (B) , and star magnitude (I) from their 

values at calibration. The cross terms, i.e. those terms involving 

in o and V in 6 , occur because the nominally orthogonal planes in 
y X 

which -jt and g are measured may be slightly nonorthogpnal, and because 
certain phenomena such as a strong axial component of the Earth's 
magnetic field cause an effective rotation of the sensor axes. 


If one considers the errors in a and B to be small compared to 
their magnitudes, it is then possible to write 


Op + Qg = (Aj^ + «A^) + {A2 + + . . . , 

= (B, + 6B,) + (B„ + + . . . . 

p e 1 1 2 2 X 
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or 


\ + 5A2V^ t 5A3Vy t . . . , 


B, 


6B3 + + . . , 


(3-11) 


Carrying the small error approximation one step further allows 
each of the coefficient errors to be written as an expansion, 
namely 


3-9 






iSAj iSi 


3A. 


9T 


6T + 


3A. 


SB 


t-t 


6B + 


3A. 


31 


51 + 


B=B. 


1=1. 


(3-12) 


Here and are the tracker temperature and star intensity at cali- 
bration ^ and B . is the axial magnetic field strength (nominally zero 

V . 

at calibration) * If one adds a first order error in the vbltage 
measurement^ terms of the form are added to the above expression 

for those coefficients that are nominally multiplied by voltages, i.e. 
A^, i > 1. An expansion similar to that in Equation 3-12 may also be 
constructed for It should perhaps be emphasized that the coeffi- 

cient errors in Equation 3-12 represent residual errors after calibra- 
tion and fitting of the measurements With selected functions. Both 
the nominal coefficients (A^, and the residual errors in the coef- 
ficients (SAj^, 6B.) are functions of the temperature, magnetic field 
and star intensity deviations from the nominal calibration values. 


3.2.6 Evaluation of Sensitivity Partial Perivatives 

To obtain the desired partial derivatives it is necessary to 
differentiate the general geometric relationship (Equation 3-7) . This 
equation is an expression of the relationship between the measurement 
vector 2 and the state x. We restate the equation as 




6 ) 


(3-13) 


The quantities on the left are each functions of the state vector x. 
Thus Tgjj is a function of the attitude, and Tq ^ and are functions 
Of the ephemeris. For manipulative purposes we may therefore express 
the left side of Equation 3-13 as a simple function of the state 


1<£) = ’'so'i'^oA . 

The measurement vector z is contained in 3^,(0, 3) of Equation 3-13, 
thus 


Sy = Sjj,[a(z) » 3 (2.) 1 . 
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If one now allows for estlmatable bias errors in temperature, star 
intensity# traeker alignment# etc,, the quantities on the right hand 
side of Equation 3-13 become functions of the expanded state vector# 
and one may write a general form for Equation 3-13 in terms of the 
state vector and measurement vector as 

f(x) = T(x) S^tatk# z ) # 0{x# z^) ] . (3-15) 

The sensitivity partial H ;=* Sz/9x can now be obtained in general form 
by implicit differentiation of Equation 3-15. Thus one obtains 

9f 3T 3^ 

— = — + T — t # (3-16) 

3x 9x 3^ 

which is further expanded to yield a solution for 3z/3x. Expanding the 
derivative in the second term on the right of Equation 3-16 gives 


3x 



3c£ 

— + 




3ct 

— + 

3^ 

33> 

> 3z 

1 — • 

3ql 

9x 

33 

3x/ 

\3a 

dz 

33 

3x> 

^ 331 


(3-17) 


Letting the terms in parentheses be abbreviated as and 
allows Equation 3-16 to be expressed as 


3f 9T 3z 

— _ — Sjp + TQ, + TQ„— , (3-18) 

3x " 3X ^ ^ ^3x 


Solving for 3^/3x produces the desired sensitivity matrix 


3z ^ . 

( 91 

1 9T \ \ 



H = — = (TQ,)“-^ 
3x ' 

' 


\ ■ 

(3-19) 


The inverted quantity TQ 2 involves only derivatives with respect to 
the measurement vector. Thus the same equation can be used to evaluate 
the sensitivity matrices for consider variables (H ) and white noise 
(A) by substituting x and v for x. For the star tracker, the same 
elements of x do not occur in more than one of the terms f, T# and S^. 
Therefore, for any x element, two of the three right-hand terms in 
equation 3-19 are zero. 


We next display some expansion steps in the evaluatibh of 
(TQ 2 ) ^ and Qj. The matrix Q 2 can be written 





!2i 

^^x 

SVy 




Wy 

383 

5Sj 




where the elements of S,j, are given by Equation 3-8. Note that the 
expanded derivatives for when written out in full, are actually 


3 u 3 $^ 33 

3 Vj 3 «jt 3 V^ 30 3 Vj 


If T is approximately orthogonal, the coefficient (TQ 2 >~'^ can be 
simplified as follows: 

(TQ2)"^ = T‘^ -5 (3-22) 

Since Q 2 is a non-square matrix it is necessary to compute the pseudo 
inverse 

■> qJ = (qI Qj)-^ qJ 


liV . C-hV . /!!3\ 

I \ 9 Zj/ \ 3 Zj/ 


2 3S- 3S, 3S- 3S» 3S_ 3S, 

— i k + — I £ + 1 

3Zj^ aZj 32 j^ azj 3Zj^ SZg 


3S^ 3Sj^ 3S2 aSj 3 S 3 383 

aZj^ SZj 8X2 ^^1 ^*2 


( es,\2 /3S,\2 /3S-. ,2 

t) * (i;!) ^ 


where is obtained by transposing Equation 3-20.. 




(3-23) 


Using Equations 3-8, 3-9, and 3-20, the first matrix in Equation 
3-23 may be re-expressed as 


2 2 2 
p + q cos a 


2 “1 

pu + qv cos 01 


. ■■■■T ■ _.T 

(Qa ^2^ 


(3-24) 


pu + qv cos a 


2 2 2 
u + V cos a 


where 


p - Aj + 2A^V^ + AgVy , 


q = + 2B^V^ + EjVy , 


(3-25) 


« = Aj + 2AjV^ + AjV^ , 


V - Bj + 2B5V^ + BjVy . 


Performing the inverse operation algebraically gives 

O O ^ 

“(pu + qv COS Qt) 

/ 

1 


T -1 
(Oj Qj) 


2 2 
cos a(pv-qu) 


2 2 2 2 
(pu + qv cos a.) P + q cos a 


(3-26) 


In terras of defined quantities. Equation 3-25, the terms of 
T 

Qj (from Equation 3-20) ate: 


3S./9V = p cos ct 


= p sin a sin B - q cos a cos B 


9S_/3V = -p sin a cos B - q cos a sin B 


9Sj^/9Vy = u cos a 


9S_/3V = u sin a cos B - v cos a cos B 

d. y 


3S_/3v = -u sin a cos 6 - v cos a sin B 

Y 


(3-27) 
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in sumtnary, the sensitivity pectials are given in general form 
by Equation 3-19. Since the variables x, x, V apply to only one of 
£, T and at a time. Equation 3-19 simplifies for any given variable 
of differentiation. For lioise. biases, etc, related to the spacecraft 
altitude we have 


-1 

H (altitude) (TQ,) ^ — . 

3x 

For star tracker alignment errors 


H (alignment) 


-(TQ2> 



(3-28) 


(3-29) 


and for tracker internal errors 


internal _ , 

H (to tracker) ■+ -Q 2 . ( 3-30) 

Since T is generally a product of a nominal rotation matrix and a 
small angle error matrix, ie. 


T 


T^T 


BT ' 


the alignment sensitivity can be written 


H (alignment) 



5.2.7 Sensitivity Example 

Suppose we wish to calculate a particular sensitivity matrix 
element. As an example we pick 3Vj^/3(6T) where ST is a temperature 
error (ST may be a bias if temperature is measured onboard the space- 
craft; or a consider type error if temperature is not measured) . Let 
6T be the first element of x or or v. Then from Equation 3-30 


D(6T) 


-ioV' Q 


I'll 
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Using Equations 3<-26 and 3<-27 one can write 


-1 




ca (pv-qu) 


2 2 2 2 
U + V CO “(pii + <jv CO ) 


2 2 2 2 
-<pu + qv CO ) p + g CO 


p CO p so S3 
-q CO o6 


u CO u so s3 
-V CO c8 


•p so c3 
-q CO sS 


>u so cB 
-V CO sB 


(3-31) 


where c and s denote cosine and sinCi respectively. Using Equations 
3-13 and 3-14, and the vector form of for a single element, gives 


<Qi>ji 




3A, 

— A 

4. 

2 

— X 


3T 

'X 

3T * 

V 

3T ^ 

V 

3T ^ 

3T y 


3A 


3T 


Wv \ 

X y/ 



3B-, 

— = + 

^®2 

— v„ + 

3B. 

— + 

®®4 2 

— iy + 

3B 2 
— -V + 

3T 

3T * 

3T y 

3T * 

3T y 


(3-32) 


Picking a = B = 0, i.e. the sensitivity at the center of the 
field of view, and multiplying out Equation 3-31 gives 

V u 0 


= 


(pv-qu) 


-q -p 0 
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Proiti Equation 3-32 


‘Oj'jX 


I'EA Terms 
Terms 
C 


therefore 


9z 


3x* 


3V. 


9 (ST) (pv-qu) 


(vZA Terms - u*ZB Terms) . 


(3-33) 


Examination of the forms of p, q, u, v, as given in Equation 3-25, 
show that normally 


*2 "" ®2 , 
Aj « Bj ^ 


A 


5 


•C < 


A 


6 



This is because the cross-coupling terms due to rotation and non- 
orthogonality are nominally small. Also Ag and Bg are small with 
AgVy << etc. These inequalities imply that 

p >> q, V >> u 


so that Equation 3-33 becomes approximately 

9V 1 • -1 9 

— c — a: _ . V • ZA Terms -(A) , 

9(dT) pv p 9T 


KEPRODUCIBaiTY OP THE 
rtWOINAL PAGE IS 
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or 


3V. 


(Jh. 

\ DT 


3A.J 3A, 3A. - 3A_ „ SA, 

3T 3T y° 3T 3T 3T 


■V V 1 

xo yo L 


X 


3 ( 6T) 


^ 2A^V^^ + AgVy^) 


T = T. 


(3-34) 

where T = is the calibration temperature, and the s'jbscript 
"o" for the voltages indicates the values assumed at a - ^ = 0. 

The behavior of Equation 3-34 can be seen more clearly by considering 
the tracker calibration to be linear without cross-coupling. Then 
Equation 3-34 reduces to 



9i6T) 


3A, 3A- 

— L + — £v 
3T 3T 


A 


2 


(3-35) 


The intuitive content of this equation is more clearly displayed by 
some manipulations. We first re-express the right side 



Substituting into Equation 3-35 and solving for <3a/3(6T> yields 


3a 3a 3Vj^ 

3(6T) 3Vjj 3(fiT) 


(3-36) 
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3.3 LANDMAHK MEASUREMENT EQUATION 


3.3.1 Introduction 

One of the primary instruments being planned for an EOS type 
satellite is the multispectral scanner. This device is a natural for 
navigation measurement using landmarks within the fi'ild of view. 

There is an important system advantage in using the device in this 
manner, and that is that landmark sightings can provide both attitude 
and ephemeris information simultaneously. For example, it is shown in 
Section 2 that scanner sightings of landmarks coupled with range and 
range rate measurements from ground stations could obviate the need 
for a star tracker onboard the spaceeraft. 

This section derives a general measurement equation for the 
landmark sightings including all of the known potential error sources. 
The equation is made deliberately general so that various error 
sources that are small are not overlooked and can be systematically 
dropped if desired when the equation is applied in practice. 

3.3.2 Basic Scanner Operation and Measurement 

The scanner moves a small field of view across the Earth’s sur- 
face at a rate designed to give full area coverage. At regular time 
intervals the sensor output is sampled to yield a sequence of effec- 
tive picture elements. If a known landmark has a distinguishable out- 
put, its position within a scan can be in principle determined by ap- 
plication of a recognition process. To this end it has been suggested 
that searchlights and mirrors be used as landmarks since they are 
nearly point sources with distinguishable spectral characteri sties - 
The searchlight and mirror types of landmarks yield a sharply defined 
signal that is essentially confined to two picture elements.' Thus if 
the approximate position of the landmark is predictable over a limited 
part of the scan » it may be possible to use a finer signal sample 

interval in order to establish the position of the landmark to within 
one optical resolution element. 
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Fondairidhtal to the navigatibh problem is the association of 
event times Within the scanner system with universal time which is 
used in ephenieris calculations. Since the scanner scans at a Known 
rate, the inertial scanner pointing direction can be obtained by db~ 
serving the universal time at the beginning of a scan on which the 
landmark is sighted, and adding to the estimated starting positiQn an 
angle which is a function of the scan rate and the time between scan 
initiation and landmark detection. The fundamental measurement is 
therefore the time interval between scan initiation and landmark 
sighting. 

If the navigation calculations are to be done at an earth facili 
ty it will be necessary to transmit from satellite to earth the scan 
initiation timing pulse (probably for several scans} and the time in- 
terval between detection and initiation. The latter can presumably be 
processed onboard the satellite. 

3.3.3 Formal Measurement and Noise 

The linearized measurement equation considering the various cate 
gories of noise is 


3Z 

3z 


&z 

6z = — 

5x + — — 

6x + 



■ 

— c 

3v 


where Xf Z defined in Equation 3-1. For the multispectrai 
scanner ^ is simply a scaler £ At. 


The geometric equations (discussed below) relating scanner lindi^ 
of-sight to orb it- inertial coordinates involve the basic measurement; 
(time of landmark sighting) by means of direction cosines of the scan 
angle a. The linear scan angle-time relationship asstimed here is only 
one of a number of possible ot(t) relationships which result from var- 
ious mechanical scanner designs. The linear relationship is relevant 
to rotating mirror scanners as in the Nimbus and ITOS satellites, the 



ERTS scanning mirror, and rotating satellites such as PiTS I,, II> 
Pioneer 10, and SMS. We write the linear (error-free) relationship as 

a = a(t ) + ctAt, (3-37) 

o 

where a(t^) is the angular deviation of the scanner field of view from 
null at scan initiation, & is the scan rate, and At is the landinark 
detection time interval. The true value of a includes errors whose 
sources are identified by examining the differential form of Equation 
3-37. Thus 


6a = 6a + (6ct) (At) + (a)6(At). (3-38) 

o 

Here 6a^ represents an tmcertainty in the scan field width, 6a repre- 
sents a scan rate uncertainty, and 6 (At) a timing error. The latter 
consists of two errors, namely 

6 (At) = 6t- - 6t (3-39) 

a o 

where t , is a time-of-detection error and 6t is a scan initialization 

(3 o 

error. Scanner related sources for these errors are listed below in 
Table 3-2. These were obtained in part from discussions with O. Wein- 
stein of GSFC and from Reference 18. The initial scan position bias 
is due partly to the internal scanner mechanism and partly to the 
alignment of the scanner with the spacecraft body axes. The scanner- 
spacecraft body alignment error is treated separately as a small angle 
rotation matrix in the discussion which follows. Angular alignment 
errors are expected to be on the order of 100 u rad. 

The consider biases are slowly changing errors of random ampli- 
tude, for example, errors associated with changes in calibrated settings 
due to temperature fluctuations. 

Since no time-of '^detection algorithm has been designed and no 
tests have been made, the random error associated with detection cannot 
be closely estimated. However the searchlight experiment data 
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TABLE 3-2 

MSS LANDMARK SKATING ERRORS 

EstimataJale Bias Errors (la) 

Inifrial Scan Position (5a ) 

o 

Scan Rate Bias (5a) 

Consider Bias Errors (la) 

Initial Scan Position Due to 

Thermal Changes and Wear (5a^) - < 4 p rad 

Scan Rate Drift (6a) - .003 % 

Initiation Time Error (5t^) - 10“® % of interya). 

White Noises ( la ) 

Side-To-Side Scan Line Jic;er (^6^) *■ 4 P rad 

Scan Line Synchronization (6t^) - 4 u ra4 

10 P rad ? 


- 4 p rad 

- O.l % 


Time of Detection (6t^) 



generally show the landmark confined to a single or at most two picture 
elements. Thus If the midpoint bf the maximum signal picture element 
is used as a marker the uneertainty will be roughly a half picture 
element. The discrimination accuracy could be improved by using a 
smaller sampling interval for a period defined by the expected landmark 
sighting time interval. 

3 , 3 , 4 Scanner- Inertia I Coordinate Relationship 

The nominal vector relationship between landmark, spacecraft and 
scanner line of sight vectors in spacecraft body Gbordinates can be 
written as 

^BO^OI— ” * (3-40J 


Here T^^ and T^^. are the rotation matrices defined above, and Tgg is a 
rotation describing the orientation of the scanner in spacecraft body 
coordinates. Previously (Ref, 1) T^g consisted of a single 90® rota- 
tion corresponding to an x axis boiesight aimed straight down along 
the body z axis. In the interest of generality T-^ is defined here to 

A Op 

be a three angle rotation matrix, u is a landmark position unit vec- 
tor in inertial coordinates and Im is a unit vector along the scanner 
line of sight towards the landmark. In the presence of errors the gen- 
eral foarm of Equation 3-40 changes only by the matrix T(6 y) which rep- 
resents small angle errors in the scanner alignment with the spacecraft 
body axes. T(5 y) is introduced as a multiplier of T--!> The error form 

OP 

of u is given by 

>, 1 - r + 

u = (3-41) 

d 

where r is the satellite position vector, ^ is the landmark position 
vector, and ^ is a landmark position error. The scalar d is just the 
magnitude of the vector numerator. Figure 3-2 illustrates the land- 
mark sighting geometry. 
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the scanher x axis £iQiht along the center of the scan 
allows the alignment matrix T^g to be defined as 



(3-42) 


The individual rotations are defined as follows; 90^ about body y axis 
through angle A about Z'; through angle B about X'; finally through an- 
gle C about X** In the scanner coordinate system at any instant, the 
scan direction will be defined by the nominal scan angle a. Thus the 

A 

nominal form of ^ is just 


'cos (a + az)> 
o 


Im = I sin [a 4* az) | , 
— ■ o 


(3-43) 


where z is the measured time of detection z = t^ - t^. Since the line 
scan type of multispectral scanner which we are considering here regis- 
ters six scan lines simultaneously, it is possible to have a displace- 
ment error in a landmark sighting when the landmark image falls between 
two detectors . Considering this displacement to be the angle S , and 
adding the error sources in a (Equation 3-38) , gives the actual ^ 
value, namely 


Im 



(a + 6a + 6az + a6Z) 
o 

(a + 5o^ + 6a2 + a6Z) 
o 

sin B 


cos B 
cos B 


(3-44) 


3-24 








The complete geometric equation in the presence Of erroirs in the mea- 
surement is thus (combining Equations 3-40, 3-41, and 3-44) 




<cos (a + Sot + Saz *»* otSz) cos S' 

'O' 


^BO^OX d ( * ^^y ~ 3^ I " ■♦■ 6a? a6?l COS 6 


1 , + 61 , - ? 

z z 


sin S 


(3-45) 


3.3.5 Evaluation of the Sensitivity Partial Derivatives 


Equation 3-45 can be written in the same general form as Equation 
3-15, namely 

f(K) = T(x)^[a(x,z) ,S(3t,?) } . 

To obtain the sensitivity to the state we again solve 


3f 3T(x)'> 31m 

— = + T(x) — 

3x 3x ” 3x 


for 3z/3x. The result is (here z is a scaler) 


H = 


3z 


3x 


= (TQ2> 


" (i - ^ ^ 


(3-46) 


Rewriting with T ^ multiplied out gives 




(3-47) 
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^call that previously 


/\ A 

( 3^ 3a 38 \ 

3a 32 36 3z / 


In the present case this becomes [since 6 = 6(z)] 



( 3 - 48 ) 


where partial derivatives are evaluated at nominal (error-free) condi- 
tions. Taking the inverse gives 



1 

— (-sin a cos a 0) . 
a 


(3-49) 


It remains to calculate the other terms in Equation 3-47. Let the 
state vector be defined as follows: 


X — ( 9 f \J) / J 


®x'®y'®z' x,y,z; 


61^,61y,61^,6a,66#iS Yf 6a) 


The angles 9, 4i, t() are the altitude definition angles contained in the 

rotation Tg^ introduced in Equation 3-40. This matrix is defined in 

Section 2.4. B , B , B are gyro bias drift rates also defined in Sec- 
X y z 

tion 2.4. The quantities 5_1 are landmark position biases corresponding 
to an uncertainty in the geoid location of a known landmark. For ex- 
ample, if the landmark is a search light surveyed in by common survey 
techhiques, the maghitude of t,his error is approximately 10 meters 
&tate vector elements x, y, z are just the satellite coordinates, Sy 
is the scanner-body alignment bias (may have three components) , and 
5a, 56 and 5a are as defined above. 


The first three state elements enter Equation 3-46 only through 
the first term in brackets. These we may write 








3T 


BO 


dx 
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Now T = Tgg. Since Tgg is orthogonal 

T~g - Tgg, and since T(6y) is infinitesimal T”^(6y) w T^(6y). Substi- 
tuting (defined in Equation 3-49) and T~^ yields 




1 

— (-sin ot cos a 
a 


^ ^OlH ' 
3x 


(3-51) 


where x here represents only elements 6# <t>, T(6y) does not appear 
because its nominal value is the identity matrix. Evaluating the par 
tial derivative gives (the nominal values of 9 and are 0°) 
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Elements of. the sensitivity vector H for the first three state elements 
are thus 


3z 


3G 


= JhgK , 


GZ 

dz 


= Jh,K 
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Where J - 


■'L s 


■'oil! 


Elements of H fOr the second set of state vector elements also 
involve only derivatives of f , thus 


-1 T 

Hjl - Q2 Tgg TgQ Tqj 


(3-53) 


Evaluating the derivatives 6£ u gives 


3u 1 

' ” hj. “ ■ 3 

3(51j^) X d 


r-„2 


-u u 
X y 


-u u 

Z X 


3u 


3(6ly) 


-u u 
X z 




-u u 

z y 


-u u 
X z 


3u 1 

^ ~ ill ” 3 ' 

3(fil_) ■‘■z d^ 


“U u 

"•y z 
2 ^ ^2 
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where u = 1 “ x etc. according to Equation 3-41. Nominally ^1 - 0. 


then the H elements associated with 61 are 


Let I ■ Q 


-1 «T 
^BS 




dz 

I h , , 

3(61 ) 


3z 

= I h, , 

3(51 ) V 


3z 

= I h, . (3-55) 

3 ( 61 ^) 


The sensitivity to satellite position is identical to the above 
except for a change of sign, thus 

3z 

— = "I h 1 * 

3x X 


3z 


— = -1 h, 
3y y 


32 

-I hi 

9(z) ^z 


(3-56) 


Sensitivity to 6y involves the second term of Equation 3-47, 


Thus 


H 



3T(6y) ^ 

Im 


3(6 Y) 


If we write 


then 


and 



The remaining three state elements involve the last term of 
Equation 3-47, and can be expressed as: 


^i(6a^,66,6a) 


-1 

-Q2^Q 


1 * 


Recall that is defined as 31m/ or in this case 


Q 


1 


( 3 1m 3tt 
3 a 3x^ 


3 ]^ 38 

m 

36 3x 



One of these terms is therefore associated with each of and 
6a. Differentiating Equation 3-44 where a is taken to be the "true" 
value, i.e. a = a + 6a, with 6a representing the error terms, gives 
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thus for this state element 



1 

— t“sin a cos a 
a 



For 63 we have 


1 

H-g = - — (-sin a cos ot 



Finally for 
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3.4 C0I4CLUSI0NS 

Xt is possible to derive general and realistic measurement equa- 
tions for optical navigation sensors following a three-step procedure. 
The steps are: 

(1) Determine the actual measured quantity in terms of the sen- 
sor operation (e.g. pulse time* voltage levels, etc.}, and 
the types of noise associated with the measurement (e.g. 
bias, white, etc.). 

(2) Write the geometric state, measurement, and error relation- 
ship in general form 

g(^) = f(x) 

where z is the measurement vector and x is the state vector 
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(3) Differentiate this relationship implicitly to obtain the 
linearized measurement equation used in the linear estima- 
tion of nonlinear systems. 


These steps have been applied to star tracker and multispectral 
scanner measurements. Specific models of these devices were considered, 
and the resulting measurement equations are fairly comprehensive. 



SECTION 4 


HORIZON UTILIZATION STUDY 


4.1 INTRODUCTION 


In this report some preliininary r^ults are 4iV^h for the Horizon 
Utilization Study, which is an investigation into the use of infra-red 
(IR) horizon measarements to estiniate satellite attitude and: gyro bias 
drift for a geosynchronous earth-observing satellite. It is assunted ^ 
;hat the satellite has three orthogonal body-mounted gyros aligned with 
it*s principal body axes> and that the satellite ihaintains d 
vertical attitude history like that used in the Modified Landmark Utili^ 
zation Study of Section 2. 


The primary objectives of this study are to review what is pre- 
sently known about the Earth's horizon in the IR region, and determine 
what might best be achieved in the estimation of satellite attitude sind 
gyro bias drift when using the most accurate instruments for horizon 
sensing. Unfortunately, there has been some difficulty irt identifying 
those instruments which might be best for this application since none 
oi the known instruments have accuracies which are comparable to the 
small angular uncertainty (approximately 0.004 degree) of the IR horizon 
when observed at geosynchronous altitudes. The typical accuracy of 
devices used for IR horizon sensing is about 0.1 degree. Design studies, 
by various companies, such as Hames Engineering ' ' ' , indicate that 

the better horizon sensors are capable of accuracies near 0.05 degree, 
with detector-preamplifier noise being one of the limiting factors* In 
addition, these instruments have usually been designed for operation at 
much lower altitudes than what is desired is this study, and must there- 
fore work with 'the much iarger angular uncertainty of the IR horizon* 
Asyet no instrument has been identified which v7ould have , a measurement 
accuracy comparable to the uncertainty of the IR horizon when observed 
at geosynchronous altitudes. However, for the purposes of this study, 
it is assumed that such an instrument does exist and has an overall 
measurement accuracy (including horizon uncertainty) of CX.Ol degree or .-, 
better. 
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In this report the performance results are given for one type of 
hypothetical horizon sensor. This sensor is assumed to be a body-fixed 
device which simultaneously senses the Earth's horizon at a number of 
equally spaced points and provides an indication of local vertical. In 
other words, it may be looked upon as a device which essentially 
tracks the center of the Earth. The measurement equations for this 
device are presented in Section 4.3.2. To simplify matters, it is 
assumed in this first analysis that the Earth's horizon is circular. 
However, as indicated later in the comments on the performance results 
(Section 4.4), the next step in this study would include consideration 
of the actual non-circular shape of the Earth's horizon, which would 
provide some additional information on yaw attitude. 

It should be noted that measurement equations have also been 
developed for a second type of horizon sensor, but are not presented 
in this report since they have not yet been incorporated in the computer 
simulation. This device consists of two small FOV sensors which spin 
about body-fixed axes in such a manner that one conically sweeps across 
the northern hemisphere while the other sweeps in a similar fashion 
across the southern hemisphere. 

4.2 INFRA-RED HORIZON UNCERTAINTY 

(22 23) 

Various studies ' have shown that the 15vi bands of CO 2 pro- 
vide horizon markers that display the least sensitivity to variations 
in atmospheric phenomena such as cloudiness and surface effects. Accord- 
ingly, most of the theoretical and experimental programs aimed at IR 

horizon fluctuation assessment have been confined to this spectral 

(24) (25) 

region. Theoretical analyses by Wark , and Thomas et al have 
shown maximum altitude fluctuations of 3 to 4 kilometers for the one- 
half maximum radiant intensity marker after corrections for latitude 
and season have been made. These fluctuations are primarily due to 
short term temperature variations. Fluctuations of this magnitude have 
been confirmed experimentally by Dodgen' ' , and Girard' ’ . For a 
satellite at 1000 km altitude, a 4 km displacement on the horizon sub- 
tends an angle of about 0 . 06 degree . At geosynchronous altitude the 
angular subtense is 0.004 degree. 
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4.3 STATE AND UEASUIffiMBNT EQUATIONS 

In this Initial covariance analysis » frequent use is made of the 
analytical techniques and definitions previously described in Section 2 
for the Modified Landmark Utilization Study. Also included in this 
study are most of the coordinate systems defined in Section 2.3. 

4.3.1 State Equations 

The state of the system x is defined as the following six dimen- 

a 

sional vector i 


where 

0 ,<|i,tp - Euler angles defining the body attitude with 

respect to the orbital frame. 

~ Bi^s drift of roll (X) , pitch (Y) , and yaw (Z) 

* y ^ gyros. 

The state equation fOr this system is the previously derived 
Equation 2'’21 of Section 2. Since this equation is non-linear, use 
is made of the corresponding linearized expression of Equation 2-28 
where the adopted state parameters represent the perturbations of the 
original parameters. The transition matrix for the linearized state 
equation is that given in Equation 2-33. 

4.3.2 Measurement Equations 

As previously indicated in Section 4.1, the horizon sensor is 
assumed to be a body-fixed device which simultaneously senses the 
Earth's horizon at a number of equally spaced points and provides an 
indication of local vertical which is analytically represented as a 
unit LGS vector to the center of the Earth. 

*For^ the purposes of this study, it is assumed that the satellite 
position and velocity are known without any uncertainty. 
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The estimated direction of this unit LOS vector in body co- 
ordinates is: 


u 


LOS 


’’bo ’^01 *-£> 


(4-2) 


where r is the estimated satellite position vector in basic inertial 


coordinates, and T._ and T__ 
OX so 

defined in Section 2.3. 


are the coordinate transformation matrices 


The unit vector defining the measured LOS to the center of the 
Earth in body coordinates is: 


- LOS 


10 0 

0 ca -sa 

0 sa ca 


c8 0 s8 

0 10 
“sS 0 cB 


(4-3) 


where a and B are Euler angle rotations about the roll (X) and pitch 
(Y) axes of the spacecraft, respectively. These measured angles con- 
sist of the true angles plus a measurement noise (v) which is assumed 
to be white and gaussian. In equation form. 


a = + V 

True a 


® ” ®True '^B 


The linearized measurement equation for this case is: 
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6z = H (Sx + A (Sv 


(4-5) 


where the matrix H is: 


H = 


LOS 


3x 


— 


-fe [fBO ■'oi <-£'] 




9x 


Tqj Unit (-r) 


l^==So 


= r Hq i i 

L • I 


0 I 0 

0 ; 0 


is] 


(4-6) 


*S~ee footnote on page 2-20. 
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(fi=i()==0) and 


where is the nominal state (i.e., 
— o 


Se 




-c6 -s0 0 

6 0 


[■ 

[- 


0 0 
c6 -s0 


Tqj Unit (-r) 


Tqj Unit (-r) 




0 0 

-30 C0 


Tqj Unit (-r) 


The matrix A in Equation 4-5 is; 


A = 


= S LOS 


3 V 


a =a 
B =6. 


0 

-cao c6^ sa^ c6^ 


which for the nominal case of a 

o 


B = 0 is; 
o 


A = 


0 1 

-1 0 


4*3,3 Filter and Smoother Equations 


(4-7) 


(4-8) 


(4-9) 


(4-10) 


(4-11) 


The computer program written for this study allows for the 
selection of either a forward filter estimation scheme or a Fraser 
Two-Filter Smoother ^ ^ . In either case, the data interval is assumed 
to begin at time tjj = 0 and ends at time t^, with measurements occur- 
ring at equally spaced times tj^, t 2 » t^ • • . r or alternatively, at 
equally spaced orbit angles • • *'^n‘ time between 

measurements is represented by Atj^ ss = constant. The time of 
interest is assumed to be t^ where tj^ £ t^ £ 


The forward filter is a standard Kalman filter using the 
following eguations to recursively update the covariance matrix of 
the state estimation uncertainties from t„ to the time of interest. 


1 

*^k , k-1 

Vi *k, k-1 * 


(4-12) 

II 

? T 

"k'’^ 


(4-13) 


!I - 

«k> '■k 'I - \ “k'^' * 

Wk 

(4-14) 


where Pj^ is the covariance of the state estimation uncertainties just 


before t,, 


is the covariance of the state estimation uncertainties 


just after tj^; Q is the covariance matrix of the gyro random drift 
(Equation 2-18) ; and G is 


G = 



(4-15) 


where is given in Equation 2-27. 


The Fraser Two-Filter Smoother makes use of two Kalman filters, 
one to process the data forward from the beginning of the data inter- 
val to a point of interest, and the other to process the data back- 
wards from the end of the data interval to the same point of interest, 
Thfe resulting estimates of the two filters at the point of interest 
are then combined in an optimal manner to obtain a smoothed estimate. 
The forward filter is the filter previously described. The backward 
filter is also a Kalman filter; however, it must first be expressed 

I 

in information form since the value of P. at end of the data inter- 

^ • 

val is unknown. To do this an information matrix is defined as 
follows : 

Uj^ = (4-16) 


where the initial value of this matrix at the end of the data interval 
is assumed to be: 


tJ 


k =n 



0 


(4-17) 
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where n corresponds to the time of the last measurement. Starting at 
t. = t_, the information matrix is updated as follows: 

Jv 'll 


Ur 


== U„ 


T 

% 




-1 


H, 


(4-18) 


Afteirwards, the matrix is propagated backwards as follows with updates 
being made at each measurement time in accordance with Equation 4-18; 


U. 


k-1 




k , kl 


[(X - G^) Ojj. (I - 


+ Jfc k-1 


(4-19) 


where 


Uj^ G (G" 


Uj^ G + 


(4-20) 


The final smoothed value ^j/n covariance matrix at the time of 

interest (tj) is obtained from the two filter values Pj and as 
follows : 


"j/n ■ « - *1 “j» ^3 - "3 "3*'' * •'3 “3 ^3 


(4-?.l) 


where 

Id - p! u!)"^]'^ (4-22) 

f I 

and Uj and P^ are given by Equations 4-19 and 4-12 respectively. 
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4 . 4 performance RESUL1TS 
4.4.1 Introduction 

The results presented in this section represent the perfori^hce 
obtained with the Fraser Two-Piltar Smoother for a satellite in both a 
geosynchronous and a sun synchronous orbit. The latter orbit was used 
primarily to provide data for comparison with the results which were 
obtained at geosynchronous altitudes. 

In all cases the satellite was assumed to be in a local vertical 
attitude with rotation occurring only in pitch. The nominal values of 
the error sources and parameters are as follows: 

Initial State Uncertainties {lo)* 

Attitude (Pitch, Roll, Yaw) - 60 arcsec (each) 

Gyro Bias Drift - 0.03 deg/hr (each) 


Horizon Sensing Errors (lo) 

Assumed to be a white noise. Values given with performance results.** 
Gyro Error Sources (1 ) 

Random Drift - 0.01 deg/hr (white noise) 

Quantization - 0.1 arcsec 

4.4.2 Performance Results for Geosynchronous Satellite 

In Table 4-1 the performance results are shown for a geosyn- 
chronous satellite. Data is presented showing the effect of variation 
in the total data processing interval (i.e., number of orbits), the 

* To simplify matters in this preliminary study, the uncertainties 
in satellite position and velocity are assumed to be zero. 

** It should be noted that the horizon errors given with the perfor- 
mance results are treated as equivalent angular errors in indi- 
cating the direction of the center of the earth. The values 
given represent the one sigma values about the pitch and roll axes. 
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FRASER SMOOO^ER PERFORMANCE FOR GEOSYNCHRONOUS SATELLITE 


Data 

Interval 

(orbits) 

Update 

Interval 

(degrees) 

Horizon 

Error 

(arcsecs) 

Estimation 

Point 

(orbits) 

State Es tiroatipn Uncertainies 


Attitude 

farcseci) 


|gQj|] 

Roll 

Yaw 

X 

■■■ 

2 

2 

20 

36 

0.25 

9.9 

10.5 

60.6 

4.4 

0.12 

0.44 

2 

2 

36 

0.25 

3.3 

3.7 

59 b 6 

4.3 

0. 05 

0,14 

2 

0.5 

36 

0.25 

1.9 

2.3 

59.5 

4.3 

0.04 

0.08 

4 

20 

36 

2 

4.6 

7.5 

59.5 

4.4 

0.05 

Qi .31 

4 

20 

2 

2 

0.8 

1.0 

59.4 

4.3 

0.02 

0.03 

8 

20 

36 

0.25 

6 . 3 

6.6 

60 

4.3 

0.02 

0.22 

8 

89 

36 

0.25 

11.9 

11.4 

61 

4.4 

0.03 

0.45 














horizon sensing error, the interval between updates with horizon meas- 
urements { in degrees of orbit) , and the point at which the smoothed 
estimate was made {i.e., number of orbits since the start of the total 
data processing interval). 

The first three rows of Table 4-1 show the effect of variation 
in update interval. Here it is seen that a decrease in the update 
interval (or an increase in the frequency and number of updates) 
improves the performance in estimating pitch, roll, pitch (Y) gyro 
bias drift, and yaw (Z) gyro bias drift. However, it is noted that 
almost no change occurs in the performance for yaw attitude and roil 
(X) gyro bias drift. 

Cases 4 and 5 of Table 4-1 show the perrormance for four orbits 
of data processing and two different values of horizon sensing error. 

It is seen that reducing the horizon sensing error from 36 to 2 arcsecs 
results in a significant improvement in all of the state parameters 
except for yaw attitude and roll gyro bias drift which remain essent- 
ially at the same values as in all other cases of Table 4-1. 

Cases 6 and 7 of Table 4-1 show the performance for eight orbits 
of data processing and two different update intervals. In comparing 
Case 6 with Case 1, which are similar except for the size of the data 
processing interval, it is seen that an increase in the data process- 
ing interval results in a definite improvement in performance except 
for the same two state parameters mentioned previously. 

The fact that no significant difference occurs in the performance 
for yaw attitude and roll gyro bias drift for all cases in Table 4-1 
does raise some question as to why this is so. At first glance, it 
would seem logical that no improvement would occur in yaw since the 
horizon measurements in the present situation provide information only 
in pitch and roll. However, it should be noted that the spacecraft is 
assumed to have a nominal local vertical attitude at all times, and 
that information on roll attitude at one point in orbit becomes infor- 
mation on yaw attitude at a point 90 degrees ahead or behind in the 
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orbit. This fact is clearly illustrated by the following two equations 
showing the relationships between the errors in roll ( 4 >) , yaw (^) , roll 
gyro bias drift and yaw gyro bias drift (Bg) :* 


642 

= C 


- (SB + 
cu 

U-CL 5B 



ID 

\H 4 J / 

6*2 

= s 

6*1 + 

6B 

- - 6b 
( 1 ) z 

(4-24) 


where the subscripts 1 and 2 denote two points in orbit separated by 
the orbit angle X; S ^ Sin (X); C - Cos (X); and ui is the orbital rate. 
If X = 90 degrees, Equation 4-24 becomes; 

S»2 = + i 6B^ - i 6B^ (4-25) 

which shows that a roll error ( 6 (I>j^) at one point in orbit becomes a 
yaw error ( 6 iJ» 2 ) 3 t a point 90 degrees ahead. It is also interesting 
to note that Equation 4-25 can be used to indicate the relationship 
between the uncertainties in yaw attitude and roll gyro bias drift of 
Table 4-1 by neglecting the small values of and so that: 

4*2 = i SB^ (4-26) 

Substituting 4.3x10 ^ deg/hr (or 4.3xl0”^ arcsec/sec) for 6 B (from 

—5 * 

Table 4-1) and w = 7.27x10 radians/sec, gives 59.1 arcsecs for 6^2 

which is close to the values given in Table 4-1. This close relation- 
ship between the performance in estimating yaw attitude and roll gyro 
bias drift has been found to be the case in most of the EOS studies 
conducted at CSDL. 

The fact that there is no change in the roll gyro bias drift per- 
formance for the various cases of Table 4-1 may also seem strange since 
one would normally expect the present type of local vertical measure- 
ments to provide more information on roll gyro bias drift than yaw 
gyro bias drift. Actually, the performance in estimating roll gyro 
bias drift is better at the beginning of data processing. This is 
clearly seen in Figure 4-1 which shows the performance of the forward 

*These equations can be obtained from the transition martix in Equation 
2-33. It should be noted in Equation 4-23 and 4-24 that the errors in 
gyro bias drift are assumed to be fixed between points 1 and 2 . 
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GYRO BIAS DRIFT RATE UNCERTAINTY { XIO”** deg/hr) 



Figure 4*1 Kaliran Filter Performance in Estimating 
Gyro Bias Drift for Ciise 2 of Table 4*1 


4 ( 


r 



filter for the first 90 degrees of orbit* The results in Figure 4-1 
are for Case 2 of Table 4-1 where an update was made after each 2 
degrees of orbital motion. Note in Figure 4-1 that the performance 
for yaw (Z) gyro bias drift eventually becomes better than that for 
roll (X) gyro bias drift. Also note that the rate of improvement in 
filter performance for yaw gyro bias drift is greater than that for 
roll gyro bias drift at the end of the data interval in Figure 4-1, 
Shown also are the smoothed estimates at 96 degrees which are the 
same values given in Table 4-1 for Case 2. 

4.4.3 Performance Results for Sun Synchronous Satellite 

For purposes of comparison with the previous results given for a 
geosynchronous satellite, some performance results are shown in Table 
4-2 for a satellite in a sun synchronous orbit like that used in the 
study of Section 2, tn contrast to the previous orbit which had an 
altitude of 35,860 km and an orbital rate of 7.27xl0"5 radians per 
second, the sun synchronous orbit has an altitude of 1,000 km and an 
orbital rate of 10"^ radians per second. In Table 4-2 it is seen that 
all of the results were generated for a data processing interval of 
four orbits. Except for the difference in orbital period (and rate) , 
the conditions used to generate the results for Case I of Table 4-2 
are the same as those for Case 4 of Table 4-1. Note that the per- 
formance in estimating pitch and roll is essentially the same for these 
two cases, while the performance in estimating yaw is much better for 
the case in Table 4-2. Also note in Table 4-2 that the general per- 
formance in estimating gyro bias drift is significantly worse than 
that of Table 4-1. This result is primarily due to the fact that 
there is much less time between successive attitude updates for the 
cases in Table 4-2 because of the higher orbital rate. 

In Table 4-2 it is seen that the performance in estimating yaw 
attitude and roll (X) gyro bias drift is essentially the same for all 
of the cases analyzed. This was ailso found to be true for all of the 
cases of Table 4-1. Also note that the relationship between the 
uncertainties in yaw attitude and roll gyro bias drift of Table 4-2 
agrees with the relationship in Equation 4-26 When the appropriate 
orbital rate to = 10“^ radians/sec is used-. 
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TABLE 4-2 

Eraser smoother performance for sun synchronous satellite 


Case 

Data 

Interval 

(orbits) 

Update 

interval 

(degrees) 

Horizon 

Error 

(arcsecs) 

Estimation 

Point 

(orbits) 

State 

Estimation Uncertainties (le) I 

Attitude 

(arcsec) 

Gyro I 
(10'- 

ias Drift 
deq/hr) 

Pitch j 

Roll 

Yaw 

X 

y 

■9HI 

1 

4 

20 

36 

2 

4.3 

7.3 

27.2 

26.9 

0.6 

4.2 

2 

4 

20 

2 

2 

0.5 

0.6 

26.8 

26.8 

0.1 

0.3 

3 

4 

10 

2 

2 

0.4 

0.5 

26.8 

26 .8 

0.1 

0.2 

4 

4 

■ 

5 

2 

2 

0.4 

0.4 

26.8 

26.8 

0.1 

0.2 

.5 

4 

2 

2 

2 

0.3 

0.4 

26 . 8 

26.8 

0.1 

0.2 
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